AD-A113  381 
UNCLASSIFIED 


INTERSPEC  INC  PHILADELPHIA  PA  F/6  20/14 

APPLICATION  OF  SELF -ORGANIZING  RANDOM  ARRAYS  TO  COMMUNICATION  S~ETC(U) 
FEB  82  F  HABER*  Y  BAR-NESS*  B  D  STEINBERG  F30602-80-C-0230 

RADC-TR-81-248  NL 


DTO  FILE  COPY  AD  A 1 1  3  3  81 


.  \ 


RADC-TR-8 1-248 

Final  Technical  Report 
February  1982 


APPLICATION  OF  SELF-ORGANIZING 
RANDOM  ARRAYS  TO  COMMUNICATION 
SATELLITES 


Interepac,  Inc. 

Prof  Fred  Haber 
Prof  Yeheskel  Bar-Ness 
Prof  Bernard  D.  Steinberg 
Paul  Chien-chang  Yeh 
Judy  Herman 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


s 


DTIC 

ELECTE 

APR  1  3  1982 


D 


ROME  AIR  DEVELOPMENT  CENTER 

Air  Force  Systems  Command 

Griff iss  Air  Force  Base,  New  York  1 3441 


82  04  13  136 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 


RADC-TR-81-248  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


WALTER  J.  BUSHUNOW 
Project  Engineer 


APPROVED:  ^ 

BRUNO  BEEK 

Acting  Technical  Director 
Communications  and  Control  Division 


FOR  THE  COMMANDER: 


Js~ 

JOHN  P.  HUSS 


Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC(DCCR)  Griffiss  AFB  NY  13441.  This  will  assist  .us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

RADC-TR-81-248  d  P  4  //  3 

J  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  Subtlttel 

APPLICATION  OF  SELF-ORGANIZING  RANDOM  ARRAYS 

TO  COMMUNICATION  SATELLITES 

5  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Technical  Report 

1  Jul  80  -  20  Feb  81 

6  PERFORMING  OTG.  REPORT  NUMBER 

N/A 

7  autmo«c#> 

Prof  Fred  Haber  Paul  Chien-chang  Yeh 

Prof  Yeheskel  Bar-Ness  Judy  Herman 

Prof  Bernard  D.  Steinberg 

8.  CONTRACT  OR  GRANT  NUMBERfsl 

F30602-80-00230 

9  PERFORMING  ORGANIZATION  NAME  AN'S  ADDRESS 

Interspec  Inc. ,  University  City-Science  Center 
3508  Market  Street 

Philadelphia  PA  19118 

10.  program  ELEMENT  PROJECT,  task 
AREA  6  WORK  UNIT  NUMBERS 

62702F 

45196311 

It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Rome  Air  Development  Center  (DCCR) 

Griffiss  AFB  NY  13441 

12.  REPORT  DATE 

February  1982 

is.  NUMBER  OF  PAGES 

123 

1  4.  MONITORING  AGENCY  NAME  a  AOORESSfif  different  from  Controlling  Office) 

Same 

IS.  SECURITY  CLASS.  ' of  thte  report) 

UNCLASSIFIED 

15a.  OECL  ASSIFICATION/  DOWNGRADING 
.  SCHEDULE 

N/A 

16.  DISTRIBUTION  statement  (of  this  Report ) 

Approved  for  public  release:  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  ebetrmet  entered  in  Block  20,  If  different  from  Report) 

Same 

IB.  SUPPLEMENTARY  notes 

RADC  Project  Engineer:  Walter  J.  Bushunow  (DCCR) 

IB.  KEY  WOROS  ( Contfnuo  on  rovorao  aide  if  neceaaery  end  Identity  by  block  number) 

Random  Arrays  Interference  Cancellation 

Satellite  Antennas  Pointing  Errors 

Spatial  Filtering  Algorithms  Interference  Cancellation  Processor 

Ring  Array 

tO.  ABSTRACT  (Continue  on  reverae  aide  ff  neceaamry  and  Identity  by  Block  number) 

-*The  investigation  of  adaptive  interference  cancelling  satellite-borne 
thinned  arrays,  capable  of  forming  narrow  beams  in  the  direction  of 
signal  sources  was  performed.  It  was  determined  that  an  antenna  con¬ 
figuration  where  the  elements  are  placed  in  the  form  of  a  ring  would  be 
most  appropriate  for  uniformity  of  the  mainbeam  pattern  and  for  good 
beamwidth  characteristics.  However,  satellite  orientation  and  ground 
location  uncertainty  limit  the  pointing  accuracy.  Therefore,  _ _ 

oo ,  :;:m„  1473 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  pace  fWhtn  Data  Eni—d) 


UNCLASSIFIED _ 

HCUWITV  CLASSIFICATION  Of  THIS  *  AQt(Whm  Dtm  Enffd) 

^utilization  of  adaption  algorithms  which  are  able  to  work  through  these 
errors  is  required.  It  was  concluded  that  by  using  this  approach,  the 
required  performance  can  be  achieved.^ 


_ UNCLASSIFIED _ 

StCU*»TV  Cl  AfttlFfC  ATlO*  O*  *“‘r  PAGEfWK#*  Omtm  Entmfd) 


PREFACE 
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Philadelphia,  PA,  under  Contract  F30602-80-C-0230, 

Job  Order  45196310,  for  Rome  Air  Development  Center, 
Griffiss  Air  Force  Base,  New  York.  The  major  portion 
of  the  work  was  performed  by  Valley  Forge  Research 
Center,  Moore  School  of  Electrical  Engineering, 

University  of  Pennsylvania,  Philadelphia,  PA,  who  was 
a  subcontractor  to  Intcrspcc  Inc.  The  report  describes 
work  accomplished  between  1  July  1980  through  20  February 
1981. 
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SECTION  I 


INTRODUCTION 

Concern  over  the  vulnerability  to  interference  of  the  uplink  of  communica¬ 
tion  satellite  relays  has  led  to  considering  adaptive  spatial  filtering  as  a 
remedy.  The  subject  of  this  study  is  the  use  of  large,  sparse,  self -organizing 
satellite-borne  arrays  for  this  purpose.  The  specific  application  to  which  this 
research  is  directed  is  a  time  division  multiple  access  (TDMA)  relay  for 
widely  dispersed  users.  Reported  here  are  the  results  of  an  initial  8-month 
study  covering  the  interval  July  1,  1980  through  February  28,  1981. 

The  adaptive  interference  cancelling  array  may  be  viewed  as  a  combination 
of  a  conventional  non-adaptive  beamforming  array  focused  toward  the  desired  signal, 
and  a  secondary  array  which  adaptively  provides  cancelling  signals  to  unwanted 
arrivals  from  other  directions.  It  is  well  established  that  such  cancelling  is 
most  effective  if  the  interferers  are  sufficiently  separated  in  angle  from  the 
desired  signal  to  put  them  on  the  sidelobes  of  the  non-adpative  beamforming 
array  pattern.  Thus,  to  effectively  withstand  nearby  interferers  requires 
that  the  non-adaptive  array  main  beam  width  be  smajl  -  that  is,  that  the 
array  be  of  large  size.  Filled  arrays  of  size  require  an  unreasonably 
large  number  of  elements  so  that  sparse  structures  are  indicated. 

With  a  sparse  array,  suppression  of  undesired  signals  will  largely 
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depend  on  the  adaptive  nulling  mechanism,  since  the  sidclobc  properties  of  such 
an  array  will  not  be  adequate.  The  number  of  elements,  which  affects  the  nulling 
power,  and  the  array  structure,  are  system  design  issues  to  be  faced.  Large 
arrays  suggest  that  array  element  locations  will  be  difficult  to  control  and 
that  means  for  accepting  element  position  variation  will  be  required.  Furthermore, 
there  is  the  question  whether  or  not  the  focusing  quality  will  be  regularly 
realizable,  taking  account  of  the  propagation  medium. 

We  further  assume  that  for  TOMA  operation  the  array  be  required  to  rapidly 
bcamswitch  among  ground  sources,  utilizing  source  location  data  supplied  to  it 
externally.  This  suggests  the  use  of  a  directionally  constrained  array  processing 
algorithm  with  capability  for  storing,  recalling,  and  rapidly  updating  element  weights. 
To  assist  in  accurate  pointing,  we  see  the  possibility  for  adapting  the  algorithm 
to  utilize  additional  available  a  priori  information  about  the  ground  sources  - 
i.c.,  information  on  signal  structure. 

These  considerations  have  led  us  to  aim  our  work  at  the  following  problems. 

(1)  What  is  a  preferred  array  configuration? 

(2)  Is  atmospheric  dispersion  apt  to  affect  beam  quality  adversely;  is 
refraction  significant? 

(3)  What  can  an  adaptive  spatial  signal  processing  algorithm  do? 

(4)  What  means  shall  be  used  to  organize  the  array  assuming  reasonable 
array  element  position  tolerance? 

At  this  stage  of  our  work,  all  of  these  questions  have  been  addressed 
with  most  of  our  results  concentrated  on  item  (3).  In  Section  2  we  give  a 
technical  review  of  this  work  with  references  made  to  detailed  technical 
reports  given  in  Appendices.  Section  3  gives  the  conclusions  of  our  present 
study  and  recommendations  for  future  work. 
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SECTION  II 


TECHNICAL  REVIEW 
a.  Array  structure 

Appendix  A  on  "The  Ring  Array,"  is  a  discussion  of  system  ideas 
conceived  for  this  application  in  the  early  stages  of  our  work.  The 
principal  conclusion  is  that  elements  deployed  with  approximately  equal 
spacing  on  a  ring  would  be  appropriate  for  its  rotational  uniformity 
of  main  beam  pattern,  and  for  its  good  bcamwidth  characteristics.  The 
ring  structure  has  the  useful  property  of  structural  symmetry,  allowing 
the  elements  to  be  fed  by  a  common  oscillator  connected  via  equal  length 
cables.  Incidental  non-uniformity  in  spacing  would  result  in  sidclobos 
typical  of  arrays  with  randomly  deployed  elements;  that  is,  the  sidclobcs 
would  be  modeled  by  a  complex  spatial  Gaussian  process  with  mean  value 
of  the  sidclobc  power  response  relative  to  that  of  the  main  lobe  equal 
to  1/N,  where  N  is  the  number  of  elements.  The  number  of  elements  will 
be  chosen  to  provide  (a)  the  requisite  system  gain,  and  (b)  the  necessary 
degrees  of  freedom  to  cope  with  a  specified  number  of  interferers. 
Practical  limits  on  system  complexity  will  also  limit  the  number  of 
elements.  If  we  assume  that  individual  elements  arc  global  coverage 
antennas  with  gain  of  20  dB,  a  value  of  20dB  +  lOlog  N  will  give  the  gain 
of  the  array.  The  average  sidclobc  response  will  not  be  too  impressive; 
lOlog  (1/N)  below  the  main  lobe.  The  system  will  have  to  rely  on  adaptive 
interference  nulling  rather  than  sidclobc  suppression.  With  N  elements 
it  is  in  principal  possible  to  cope  with  (N-i)  interferers.  It  should  be 
pointed  out  that  a  large  linear  filled  array  would  require  an  unreasonable 
number  of  elements. 
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If  Che  method  of  self-organization  described  in  Appendix  G  is  used, 
a  conventional  antenna  of  about  1°  bcamwidth  (45dB  gain)  centered  on  the 
circular  array  structure  will  be  suitable  as  a  reference  source.  It  would 
be  advisable  to  provide  for  simultaneous  mechanical  steering  of  the  reference 
antenna  and  array  toward  the  ground  beacon.  This  is  suggested  in  order 
to  minimize  the  range  of  delay  differences  across  the  array,  hence  to 
maximize  the  array  bandwidth. 


b.  Medium  characteristics 

At  this  stage  of  our  work  a  study  on  the  subject  of  medium  effects 
has  been  initiated  with  a  literature  research  and  a  review  of  the  most 
pertinent  publications  (1,  2,  3,  4,  5J.  We  arc  concerned  with  (1)  refraction 
effects  which  may  bend  the  beam  as  it  enters  the  troposphere,  thus  complicating 
the  steering,  and  (2)  dispersion  effects  which  will  dcfocus  the  beam  and 
possibly  affect  the  array  adaption  process.  The  first  of  these  is  not 
viewed  as  a  serious  problem  unless  the  refraction  effect  is  severe.  Satellite 
orientation,  ground  location  uncertainty  and  other  vagaries  arc  apt  to 
limit  point  accuracy. 
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We  are  therefore  looking  to  adaption  algorithms  which  will  be  able  to  work 
through  these  errors,  as  discussed  in  2(c)  below.  Dispersion  is  potentially 
more  serious  if  it  is,  in  fact,  found  to  occur  to  excess.  The  conventional 
adaption  processes  all  are  based  on  reception  from  point  sources.  It  is 
known  that  highly  correlated  multiple  arrivals  from  different  directions 
confuse  the  array  when  it  functions  as  a  direction-finder  [6].  It  may 
pose  similar  problems  as  a  receiver,  and  the  question  of  how  it  will  behave 
under  such  circumstances  is  proposed  later  for  continued  study. 

A  recent  paper  by  Cox  and  Arnold  [1]  reports  the  results  of  propaga¬ 
tion  experiments  using  emissions  from  the  COMSTAR  satellite  at  19  and 
28  GHz.  For  our  purposes  it  is  important  to  know  whether  or  not  the  phase 
front  arriving  at  the  array  is  planar  over  the  extent  of  the  array.  The 
experiment  reported  measured  instead  the  amplitude  and  phase  as  a  function 
of  frequency  using  a  7  meter  dish  as  receiving  antenna.  Spatial  and 
spectral  measurements  in  a  multipath  medium  are  related,  so  that  such  an 
er  .iment  should  provide  clues  to  the  spatial  integrity  of  the  wavefront. 
However,  making  the  measurements  with  a  large  dish  has  the  effect  of 
averaging  over  the  wavefront;  amplitudes  at  different  frequencies  would, 
as  a  result,  not  show  great  variation,  though  phase  shift  measurements  should 
show  a  tendency  to  randomness  with  frequency  if  the  phasefront  Is  distorted. 

The  results  reported,  however,  that  no  significant  phase  distortion  was  observed 
over  a  band  of  500  MHz,  suggest  no  corresponding  wavefront  distortion  over 
the  extent  of  the  antenna.  Difficulties  with  the  phase  measurement, 
conceded  by  the  authors,  leave  questions  about  this  conclusion.  Also, 
other  papers,  such  as  that  by  Harris  and  Hyde  [2],  though  not  explicit 
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about  the  mechanisms,  lead  one  to  believe  that  phenomena,  other  than  merely 
rain  attenuation,  play  roles.  Because  of  the  importance  of  this  matter 
to  the  functioning  of  the  adaption  process,  we  believe  it  essential  to 
carry  out  a  definitive  experiment  to  determine  the  spatial  integrity  of 
the  wavefront,  preferably  by  direct  simultaneous  measurement  at  a  number 
of  spatially  distributed  points. 

c .  The  Spatial  Filtering  Algorithm 

The  preliminary  candidates  for  the  processing  algorithm  were 
(1)  the  Widrow-Compson  scheme  usivsg  an  internally  generated  reference 
[8,  9],  and  (2)  the  Applebaum  directionally-constrained  scheme  f7).  The 
former  depends  on  a  knowledge  of  some  unique  signal  property;  in  this  case, 
the  signal  is  assumed  sent  using  a  spectrum  spreading  code  known  to  the 
receiver.  The  latter  depends  on  a  knowledge  of  the  direction  of  arrival 
of  the  wanted  signal.  In  our  application,  the  directionally-constrained 
array  appeared  to  be  a  natural  choice  since  knowledge  of  source  location, 
perhaps  to  within  some  tolerance,  is  known.  As  a  first  step,  we  developed 
an  analysis  of  a  version  of  the  Widrow-Compton  scheme  (using  a  reference¬ 
generating  loop  without  limiter  and  with  one  element  unweighted),  as  dis¬ 
cussed  in  detail  in  Appendix  B,  and  a  comparison  analysis  of  this  scheme 
with  the  Applebaum  directionally-constrained  approach,  as  discussed  in 
detail  in  Appendix  C.  Both  analyses  were  done  for  a  2-element  array  in 
order  to  obtain  the  properties  without  undue  complication.  Briefly,  the 
comparison  showed  that  under  ideal  conditions  (reference  perfectly  produced 
in  the  first  and  zero  pointing  error  in  the  second),  the  first  does  a 
better  job  against  interference  close  to  the  desired  signal,  and  the  latter 
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does  a  better  job  away  from  the  desired  signal.  The  effects  of  potential 
errors  in  either  the  reference  or  the  pointing  direction  remained  to  be 
accounted  for.  The  particular  concerns  are  (1)  uncompensated  phase  shift 
in  the  reference  loop  in  the  first  scheme,  and  (2)  simply  not  pointing 
'accurately  in  the  second  scheme. 

In  connection  with  the  effect  of  pointing  error  in  the  directionally- 
constrained  array,  we  carried  out  calculations  of  Signal-to-Interference 
plus  noise  ratio  (SINR)  vs.  pointing  error  for  a  variety  of  conditions: 
numbers  of  elements,  numbers  of  interferers ,  and  interference  arrival 
directions.  The  sensitivity  to  pointing  error  turned  out  to  be  substan¬ 
tial;  in  the  worst  case, 0.1°  pointing  error  resulted  in  about  3  dB  loss 
of  SINR.  The  detailed  results  are  presented  below  in  Appendix  D. 

This  work  was  carried  out  with  one  realization  of  a  linear  array 
with  arbitrarily  placed  elements.  We  digressed  to  carry  out  a  simulation 
with  various  array  element  layouts  randomly  chosen  to  determine  the  effect 
of  layout  variation  on  the  SINR  with  zero  pointing  error.  These  results 
are  presented  below  in  Appendix  E. 

Generally,  for  interference  away  from  the  main  beam,  the  SINR  (for  the 
7-element  array  simulated)  varied  about  3  dB,  though  we  believe  that  the 
variation  will  be  less  for  arrays  with  larger  numbers  of  elements. 

The  conclusion  that  pointing  error  of  relatively  small  magni¬ 
tude  causes  significant  loss  of  signal,  suggests  the  need  for 
methods  which  overcome  the  loss.  In  Appendix  F,  a  brief  report  is  pre¬ 
sented  of  our  first  attempts  in  this  direction.  Two  alternatives  are 
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proposed.  The  first  is  a  hybrid  of  the  Applebaum  and  Widrow-Compton 
approaches.  The  second  approach  adds  phase  compensation  to  the  reference 
loop  and,  more  to  the  point,  corrects  the  steering  vector.  Since  the 
preparation  of  this  report,  additional  analysis  has  been  carried  out  on 
the  hybrid  (the  first  approach)  which  does  in  fact  show  that  the  sensitivity 
to  pointing  error  is  diminished.  Qualitatively,  we  view  the  operation  of 
this  scheme  as  follows:  The  approximate  steering  vector  puts  the  quiescent 
beam  close  to  desired  signal,  giving  a  fair  SINR  quickly.  The  reference 
generator  loop  then  has  an  adequate  signal  to  function  rapidly  to  raise  the 
SINR  to  optimum.  Though  Appendix  F  shows  two-element  arrays,  we  see  no 
problem  in  extending  the  scheme  to  multi-element  arrays.  In  particular, 
one  can  imagine  the  constrained  array  shown  in  Figure  2  of  [7]  with  the 
output  going  to  the  input  of  the  reference  loop  (as  shown  in  Figure  1  of 
Appendix  F) ,  and  the  output  error  signal  of  the  reference  loop  used  as  the 
residue  feedback  in  Figure  2  of  [7], 

Work  on  compensation  of  the  reference  loop  for  incidental  phase 
shift  has  been  carried  out  in  connection  with  another  study  done  here  [10]. 
This  work,  done  for  a  loop  with  limiter,  is  now  being  continued  for  a  loop 
without  limiter. 

d.  Array  organization 

If  one  is  to  use  a  directional  constraint  to  form  a  quiescent  beam, 
as  discussed  in  2(c)  above,  one  requires  (1)  either  element  position 
information,  or  (2)  a  reference  beacon  on  which  to  focus  the  array.  Both 
of  these  approaches  have  long  been  subjects  of  analysis  at  this  laboratory 
in  connection  with  other  studies.  In  this  application,  it  may  be  advisable 
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to  incorporate  both  schemes  as  protection  against  failure  of  either  one. 
Self-contained  systems  based  on  distance-measuring  infrared  lasers  appear 
to  be  feasible.  We  have  not  pursued  the  details  of  such  a  scheme  during 
this  phase  of  our  study.  We  have,  however,  given  consideration  to  the  tech¬ 
nique  of  focusing  the  array  on  a  beacon  and  steering  the  beam  so  formed 
toward  the  desired  source.  The  problem  we  see  is  that  the  initial  beam¬ 
forming,  if  done  by  conventional  methods  111]  in  the  face  of  strong  inter¬ 
ference,  will  be  imperfect.  This  suggests  that  the  initial  beamforming 
should  itself  be  done  using  an  interference-cancelling  algorithm.  An 
approach  of  this  kind  with  its  theoretical  basis  is  given  in  Appendix  G 
below.  It  is  based  on  a  reference  signal  obtained  from  the  beacon  (which 
is  assumed  located  far  from  interfering  sources),  through  a  separate 
conventional  antenna,  focusing  the  array  using  the  Widrow  algorithm. 

The  element  weights  developed  contain  the  steering  information  and  this 
is  isolated  by  appropriate  processing.  Once  this  is  known,  the  beam  may 
be  further  steered  toward  the  desired  source.  The  concept  remains  to  be 
analyzed  for  its  sensitivity  to  measurement  errors  and  for  the  complexity 
which  it  will  add. 

We  point  out  that  the  use  of  a  separate  reference  link  to  the  beacon 
may  also  be  useful  for  the  purpose  of  transmitting  spread  spectrum  codes, 
timing  signals,  control  information,  etc.,  to  the  satellite. 


3.  CONCLUSIONS  AND  RECOMMENDATIONS 

During  this  initial  study,  we  have  examined  the  large,  sparse,  self- 
organizing  array  for  its  potential  as  a  spatial  filter  to  be  used  on  a 
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communication  satellite  uplink.  We  have  focused  on  several  problem 

areas  -  array  structure,  effect  of  propagation  medium,  processing  algorithms, 
and  array  self -organization.  As  a  consequence  of  this  work,  we  conclude  that 

(a)  a  system  with  a  ring  structure  of  small  earth  coverage  elements 
centered  on  reference  antenna  of  conventional  form, 

(b)  a  Widrow  interference-resistant  processor  for  beamforming  on  a 
beacon,  and 

(c)  a  hybrid  interference-cancelling  processor  (see  2(c))  tolerant  to 
pointing  error  for  focusing  on  the  desired  source, 

can  be  expected  to  serve  the  intended  purpose. 

Our  work  has  proceeded  to  a  point  of  reasonable  confidence  in  the  approach, 
though  certain  problems,  as  enumerated  below,  require  attention.  Accordingly, 
we  recommend  the  following  steps  be  carried  out  in  a  continuing  study. 

(a)  Obtain  definitive  Information  on  the  behavior  of  the  propagation 
medium,  preferably  by  direct  experiment  using  spatially-distributed  sensors 
receiving  from  a  spaceborne  point  source,  and  translate  the  result  into  effect 
on  the  adaption  process. 

(b)  Continue  to  pursue  the  hybrid  and  the  self -correcting  interference¬ 
cancelling  processors  described  here,  to  establish  their  SINR  properties  and 
their  convergence  characteristics.  Develop  results  for  multi-element, 
multi-interference  cases. 

(c)  Analyze  the  method  described  here  of  beamforming  on  a  beacon  in 
the  presence  of  interference,  assessing  its  complexity  and  determining  the 
accuracy  with  which  it  generates  the  steering  vector. 
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(d)  Continue  to  study  the  effects  of  a  non-ideal  reference-generating 
loop,  and  seek  means  for  overcoming  these  effects. 

(e)  Study  the  stability  problems  of  the  various  feedback  loops  used. 

(f)  Investigate  the  problems  arising  from  the  multiple  access  feature 
of  the  system.  Consider  the  alternatives  of  beam  switching  among  users 
and  the  use  of  continuous  multiple  directional  constraints  which  allows 
simultaneous  beams  to  be  formed  on  several  users. 

(g)  Examine  the  alternative  spread  spectrum  techniques  -  frequency 
hopping  and  PN  code  modulation  for  their  respective  advantages  in  the  systems 
proposed . 
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THE  RING  ARRAY 


Let  the  largest  dimension  of  the  two-dimensional  aperture  be  called  L, 
let  the  wave-length  be  X,  and  let  the  following  three  assumptions  be  made: 
First,  the  shape  of  the  radiation  pattern  must  be  indeoendent  of  steering 
direction  throughout  the  field  of  view  that  encompasses  the  earth.  Second, 
jamming  protection  from  as  close  to  the  beam  axis  as  possible  must  be  obtained. 
Third,  the  array  is  to  be  sufficiently  large  so  that  rigidity  in  its  structure 
may  not  be  assumed. 

A  concomitant  of  the  first  assumption  is  a  requirement  for  rotational 
symmetry  of  the  array.  The  second  assumption  places  two  requirements  upon 
array  design,  one  with  respect  to  the  cross-section  of  the  main  lobe  and  one 
with  respect  to  the  side  radiation  pattern — since  these  two  requirements  are 
uniquely  associated  with  the  second  assumption.  The  first  requirement  is  that 
the  main  lobe  have  as  small  a  beamwidth  as  possible  while  the  second  require¬ 
ment  is  that  all  side  lobes  be  adequately  low  to  reject  jamming  or  that 
interference  cancellation  techniques  be  used  against  the  jammers,  or  both. 

A  circular  array  with  a  central  reference  element  is  an  attractive 
candidate.  This  array  satisfies  the  circular  symmetry  requirement  while  also 
providing  the  smallest  beam  cross  section  consistent  with  a  symmetrical  design. 
While  the  side-lobe  pattern  of  a  circular  array  normally  is  not  particularly 
attractive  relative  to  what  is  obtainable  from  a  carefully  designed,  tap.ered 
aperture,  the  third  assumption  (of  nonrigidity)  precludes  such  control  of 
sidelobes  irrespective  of  aperture  design.  Thus,  the  side  radiation  properties 
may  degenerate  into  those  of  the  random  array  regardless  of  the  designers' 
wishes.  Since  the  random  array  sidelobe  properties  are  inadequate  for  this 
purpose  (average  sidelobe  cover  level  is  1/N  relative  to  main  lobe),  inter¬ 
ference  cancellation  techniques  must  be  applied.  In  addition,  because  of  the 
nonrigidity  assumption,  some  form  of  self-cohering  using  one  or  more  pilot 
signals  from  earth  is  proposed.  The  central  element  in  the  array  will  be  the 
reference  element  both  for  adaptive  beamforming  and  for  interference  cancel¬ 
lation  , 
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The  beamwidth  of  a  continuous  ring  aperture  is  approximately  0.71X/L. 

By  way  of  comparison  the  beamwidth  of  a  uniformly  excited  rectangular  aperture 
of  length  L  is  0.88X/L.  The  beamwidth  of  a  triangularly  weighted  aperture  of 
the  same  length  is  1.27X/L. 

The  ring  array  appears  to  have  an  additional  virtue  with  resnect  to  the 
interference  cancellation  function  of  the  system.  Interference  cancellation 
is  a  nulling  process  whereby  a  null  is  formed  in  the  radiation  pattern  in  the 
direction  of  each  jammer.  The  width  of  the  null,  or  the  local  extent  of  the 
effect  upon  the  radiation  pattern  of  the  null,  is  determined  by  the  separation 
of  the  antenna  elements  from  which  the  signals  are  derived  that  form  the  null. 
By  placing  all  elements  on  the  circumference  of  a  circle  the  maximum  number  of 
pairs  of  elements  having  the  maximum  seoaration  is  obtained. 

Since  the  array  is  assumed  to  be  nonrigid  the  statistics  of  the  side 
radiation  pattern  will  be  independent  of  the  exact  locations  of  the  elements, 
provided  that  the  element  density  is  approximately  uniform  around  the  ring. 
Nevertheless,  the  designer  has  the  choice  of  equal  angular  spacing  between 
elements  or  some  statistical  distribution  about  equal  angular  spacing.  The 
former  appears  more  desirable  from  the  point  of  view  of  practicality:  given 
a  random  distribution  of  locations  the  signal  processor  must  compute  phase 
shifts  for  beam  steering  for  each  of  the  elements.  On  the  other  hand,  with 
equal  angular  spacing  between  elements  a  much  smaller  set  of  calculations  is 
required,  making  use  of  angular  symmetry  within  the  array.  The  value  of 
angular  symmetry  is  weakened,  however,  by  the  nonrigid  frame  of  the  array; 
nevertheless,  the  array  at  worst  will  be  only  somewhat  distorted,  relative 
positions  of  adjacent  elements  being  at  most  slightly  perturbed,  and  therefore 
the  underlying  angular  symmetry  may  still  be  useful. 

Two  generally  different  approaches  may  be  considered  for  the  use  of  the 
element  signals  for  interference  cancellation.  The  first  is  a  fully  adaptive 
mode  in  which  the  array  output  (or  an  error  signal  derived  from  it)  is 
correlated  with  each  of  the  element  signals,  and  the  correlation  products 
applied  as  gain  controls  to  the  element  weights.  This  fully  automatic  technique 
has  the  advantage  of  distributing  the  available  nulls  of  the  array  in  their 
optimum  locations  so  as  to  minimize  the  overall  interference  field.  In  one 
simple  scheme  successfully  used  in  radar,  called  coherent  sidelobe  suppression, 
the  signal  from  a  single  element  is  (complex)  gain-controlled  and  added  to  the 
array  output  to  remove  the  interference  from  a  single  jammer.  However,  the 


technique  has  certain  limitations  when  the  desired  signals  and  the  interference 
signals  occupy  a  wide  dynamic  range  and  when  the  interferers  are  amplitude 
modulated.  This  suggests  that  an  alternative  approach  involving  human  interac¬ 
tion  from  the  ground  via  a  data  link  may  be  appropriate.  Both  approaches  bear 
examination.  Thus,  it  is  possible  that  resilient  automatic  techniques  can  be 
designed  for  element  selection,  based  upon  array,  earth  station,  and  earth 
jammer  geometry,  as  well  as  on  automatic  complex  gain  control.  It  is  also 
possible  that  operator  interaction  in  the  selection  of  the  element  for  the  jammer 
suppression  and/or  the  adjustment  of  the  complex  weight  might  be  more  satisfactory 
under  conditions  of  severe  jamming.  A  fully  automatic  system  with  an  opera¬ 
tor  override  capability  is  likely  to  be  the  most  desirable. 

The  nonrigidity  of  the  array  requires  that  it  be  self-adaptively 
organized.  One  likely  means  is  to  self-cohere  it  retrodirectively  upon  the 
signal  source  on  the  ground.  The  same  techniques  can  be  applied  with  a 
separate  set  of  phase  shifters  to  a  jamming  source.  The  output  of  the  second 
beam  is  that  of  the  high  gain  array  pointed  at  the  jammer.  If  this  signal  is 
now  added  through  a  complex  weight  to  the  main  array  output,  the  result  multi¬ 
plied  with  the  jammer  signal  beam  output  and  sent  back  via  LMS  control  to  the 
complex  weight,  the  system  output  becomes  deprived  of  the  jamming  signal.  In 
effect  the  system  radiation  pattern  is  that  of  the  adaptively  formed  beam  upon 
the  signal  source  plus  the  complex  weighted  adaptively  formed  beam  upon  the 
jammer,  where  the  complex  weight  causes  the  sum  to  be  zero  in  the  direction  of 
the  jammer.  This  technique  uses  all  the  resources  of  the  array  both  for  main 
beam  formation  for  the  signal  as  well  as  for  jammer  suppression. 

In  summary,  one  contemplated  array  design  is  a  ring  of  diameter  L  of  N 
uniformly  spaced  elements,  plus  a  central  reference  element  used  for  adaptive 
beamforming  and  nulling.  The  beamwidth  is  approximately  A/0.7L.  A  pilot 
signal  from  the  ground  station  will  phase  synchronize  the  array  so  as  to  form 
a  retrodirective  beam  on  the  source.  Automatic  or  semi-automatic  interference 
cancellation  circuits  will  suppress  jammers. 
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INTRODUCTION 


The  adaptive  array  processing  has  become  a  way  of  life  in  the  design 
of  advanced  communication  systems,  particularly  for  automatic  beam  steering 
and  interference  cancellation.  The  Widrov  LMS  [1]  algorithm  is  widely  used. 
To  obtain  the  "reference  signal"  required  for  the  realization  of  this 
algorithm,  Compton  [2]  suggested  processing  the  array  output  using  what  is 
called  a  reference  signal  loop  (Figure  1) .  This  processing  is  intended  to 
extract  a  good  version  of  the  desired  signal  while  altering  or  suppressing 
the  interference.  For  this  to  be  possible  the  desired  signal  and  the  inter¬ 
ference  must  differ  in  some  way.  For  example,  if  the  signal  is  of  the 
SDread  spectrum  type,  then  multiplying  the  array  output  by  the  right  code 
will  spread  the  interference  bandwidth,  collapse  the  desired  signal  spectrum 
and  reveal  the  condition  of  narrowband  signal  and  broadband  interference. 
This  is  obviously  the  case  if  the  local  code  generator  is  in  synchronization 
with  the  code  of  the  incoming  desired  signal. 

However,  code  acquisition  can  not  be  obtained  when  the  interference-to- 
desired  signal  ratio  is  large.  In  this  case  the  reference  signal  loop  is 
ineffective  (produces  zero  output)  and  the  array  might  shut  itself  off  by 

[1]  B.  r-Tidrow,  J.  McCool  and  J.  Ball,  "The  Comdex  IMS  Algorithm,”  Proc. 
IEEE,  Vol.  63,  April  1975,  pp .  719-720. 

[2]  R.T.  Compton,  Jr.,  "An  Adaptive  Array  in  Soread-Spectrum  Communication, 
Proc.  IEEE,  Vol.  66,  March  1978,  p.  289. 
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turning  Che  weights  down  to  zero.  Therefore,  for  such  a  condition  of 
interference-to-desired  signal  ratio,  it  might  be  preferable  to  use  one  of 
the  power  inversion  schemes  of  interference  cancelling  such  as,  for  example, 
the  Widrow-based  power-inversion  scheme  (one  array  element  unweighted).  As 
a  result  (even  though  the  reference  loop  is  initially  ineffective)  the 
interference  will  be  suppressed  below  the  desired  signal  level  (power  inversion) , 
enabling  code  synchronization  and  finally  the  effective  reference  signal 
loop  ODeracion.  Making  the  array  control  loop  slower  than  the  code  acquisition 
time  will  assure  that  the  array  will  not  suppress  the  desired  signal  when  its 
power  at  the  input  is  larger  than  the  interference. 

Based  on  this  argument  it  is  plausible  to  assume  that  the  array 
arrangement  of  a  Widrow-based  power  inversion  scheme  with  reference  signal 
loop  will  better  handle  the  ease  of  spread  spectrum  desired  signal  than  would 
the  regular  LMS  arrangement,  particularly  when  both  the  interference  and  the 
signal  have  large  dynamic  range.  It  might  be  important  to  emphasize  that 
this  array  arrangement  does  not  require  the  limiter  which  is  very  crucial  in 
Compton's  arrangement  of  Figure  1. 

The  purpose  of  this  work  is  to  analyze  this  new  interference  canceler 
arrangement  and  consider  its  performance  with  different  system  conditions; 
namely,  when  the  system  is  in  the  code  acquisition  or  the  tracking  mode.  This 
will  be  done  by  using  the  different  bandwidth  relations  of  the  desired  signal 
and  interference  in  comparision  to  the  reference  loop  filter  bandwidth 
obtained  when  the  local  spread  spectrum  code  is  or  is  not  in  synchronism  with 
the  incoming  code..  Obviously,  when  the  system  is  in  the  acquisition  mode 
(the  local  code  is  not  in  synchronism),  both  desired  and  interference 
signals  will  be  spread  in  bandwidth  and  will  be  rejected  by  the  comparably 
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smaller  reference  signal  loop  filter  bandwidth,  making  this  loop  ineffective. 

The  desired  signal,  however,  will  benefit  from  the  resulting  power  inversion 
scheme,  pushing  the  system  toward  synchronism.  On  the  other  hand,  when 
the  system  is  in  the  tracking  mode ,  the  desired  signal  bandwidth  becomes 
smaller — and  the  interference  signal  bandwidth  larger — than  the  reference 
signal  loop  filter  bandwidth,  so  that  a  moderately  clean  replica  of  the 
desired  signal  will  be  extracted.  In  this  mode  the  array  arrangement 
operates  as  a  regular  LMS  array.  Rather  chan  assuming  only  an  ideal  filter 
in  the  reference  signal  loop,  in  our  analysis  we  will  utilize  the  filter 
responses  to  the  desired  and  interference  signals.  Depending  on  the  relative 
bandwidth  conditions,  we  will  apply  reasonable  practical  approximations  to 
these  responses  and  obtain  the  effect  of  filter  design  on  the  array  performance. 

SYSTEM  REPRESENTATION  AND  ANALYSIS 

The  array  arrangement  is  depicted  in  Figure  2.  Here,  for  the  sake  of 
simplicity,  we  use  a  two  element  array  (one  unweighted).  Generalization  to 
a  multielement  case  seems  easy.  The  signals  at  the  two  elements  are  summed 
after  being  weighted  by  a  complex  gain,  W.  The  reference-signal  loop 
comt>osed  of  a  multiplier  (spread-spectrum  decoder),  a  narrowband  bandpass 
filter  and  a  second  multiplier  (encoder).  The  purpose  of  the  first  multiplier 
is  to  decode  the  spread-spectrum  signal  at  the  output  of  the  array,  resulting 
in  a  narrowband  desired  signal  (bandwidth  of  the  data) ,  and  a  bandwidth 
spread  interference.  In  this  case  the  loop  filter  will  sufficiently  suppress 
the  interference,  leaving  an  output  which  is  a  sufficiently  clean  replica  of 
the  desired  signal.  Since  the  reference  must  carry  the  same  code  as  the 
desired  input  signal,  a  second  multiplier  recoding  the  information  signal 
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is  needed.  Notice  that  as  a  result  of  having  one  element  unweighted  there 
is  no  need,  in  this  array  processor,  for  a  limiter  in  the  reference  signal 
loop  as  was  needed  in  Compton's  processor. 

To  simplify  the  analysis  we  will  consider  two  modes  of  processor 
operation.  In  the  first  we  assume  that  throughout  the  reference  loop  the 
interference  signal  bandwidth  is  very  wide  compared  to  the  desired  signal 
bandwidth  (data  bandwidth) .  This  corresponds  to  the  case  when  the.  local 
spread  spectrum  code  is  in  synchronism  with  the  incoming  code.  Obviously, 
we  will  choose  the  filter  of  the  reference  signal  loop  to  be  wider  than 
the  desired  signal  bandwidth,  but  narrower  than  the  interference  signal 
bandwidth.  In  the  second  mode  we  will  assume  that  the  bandwidths  of  both 
the  desired  and  the  interference  signals  are  much  wider  than  that  of  the 
reference  loop  filter  bandwidth.  This  corresponds  to  the  case  when  the 
local  spread-spectrum  code  is  not  in  synchronism  and  the  code  synchronization 
system  is  in  the  acquisition  mode. 

Let  the  complex  envelope  of  the  input  signals  at  the  main  (unweighted) 
elements  of  the  array  be  designated  by  XqU)  and  S^(t)  respectively. 

Also  let 


Xt(t)  -  Ii(t)+Si(t)+Ni(t) ,  2  -  1,2,  (1) 

where  the  interference  signal  I  and  the  desired  signal  S  are  uncorrelated 
plane  waveform  signals  (If  they  are  CW  signals  then  we  assume  that  their 
frequency  difference  is  larger  than  the  array  weight  control  loop  bandwidth) . 
The  noises  processes  are  white  and  independent.  That  is. 


j8  T 

I£(t)  -  | It(t) ! e 

J9, 

S4(t)  -  |si(t)|e 
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where  9.  -9C  -  0  -  -  sintji  ,  d  is  the  distance  between  the  array 

1  b0  S10  c  S 

elements,  is  the  radian  frequency  of  the  desired  signal,  c  is  the  speed 

of  electromagnetic  waves,  and  i|/g  is  the  direction  of  arrival  of  this  signal 

with  respect  to  the  array  broadside  direction.  Similar  definitions  apply 

for  9_  .  Obviously,  depending  on  the  signal's  information  contents,  0, 

M  H 

and  0_  are  functions  of  time.  Only  when  they  are  CW  signals  they  turn  out 
bi 

to  be  constants. 

In  complex  notation  the  output  of  the  array 

V(t)  -  X1(t)W1(t)  +  XQ(t)  (3) 

The  output  of  the  reference  loop  is  termed  Xr(t)  and  the  error  signal 
becomes 

e(t)  -  V ( t)  -  Xr(t) 

-  X1(t)W1(t)  +  XQ(t)  -  xr(t) 

—  2 

For  the  LMS  processing  |e(t)[  is  minimized  by  changing  the  complex  weight 
W^(t).  Using  the  steepest  descent  algorithm  this  leads  to  the  Widrow- 
Hopf  algorithm  for  controlling  W^(t).  That  is, 

dW1(t)/dt  -  KX*(t)[X1(t)W1(t)+X0(t)-Xr(t)J  (4) 

where  K  is  the  system  gain  and  *  stands  for  complex  conjugate.  In  a 
practical  system  we  may  realize  this  by  using  a  low  pass  filter  Instead 
of  the  ideal  integration  required  by  (4).  That  is,  the  weight  is  determined 
from 

IdW1(t)/dt+[l+K,'X1(t>!2JW1(t)  -  -KX*  (t)[X0(t)-Xr(t)  J  (5) 

In  the  steady  state  (5)  gives 
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Wl(t)  -  -KX1*(t)[X0(t)-Xr(t)]^l+K|X1(t)|2]  (6) 

where  the  bar  stands  for  the  time  average  of  the  corresponding  terms. 

That  is,  for  example,  using  (1)  and  (2)  we  get 

|xi(t)  1  2  -  ll^2  +  (SjJ2  +  oN2  (7) 

X*(t)(X0(t)-Xr(t))  -  IQI  *  -IrI*  +SQS*  -SrS  *  (8) 


where  we  define  Xf(t)  -  Ir(t)+Sr (t)+Nr (t) ,  the  subscript  r  signifying  the 

reference  loop  output,  and  where  we  used  the  uncorrelation  property  between 

pairs  of  interference  signal,  desired  signal,  and  noise.  We  also  assumed 

that  the  noise  at  the  output  of  the  reference  loop  Nf(t)  is  uncorrelated 

with  N^(t) .  Notice  that  in  (7)  and  (8)  in  our  notation  we  drop  the  dependence 

on  t.  We  will  also  do  so  in  the  sequel  whenever  it  is  clear.  However, 

even  though  1^  and  S ^  are  time  dependent  due  to  their  zero  bandwidth,  the 
*  * 

terms  IqI^  and  soSi  are  time  independent.  They  are  complex  numbers  whose 

phases  are  constant  depending  on  the  direction  of  arrival  of  these  signals. 

This  is  true  if  the  distance  between  elements  is  small  enough  for  the 

given  rate  of  information  these  signals  bear.  Similarly,  we  will  also 

assume  that  the  reference  signal  loop  group  delay  is  sufficiently  small 
*  * 


so  that  and  SrS^  are  constant  complex  numbers. 

Using  (7)  and  (8)  in  (6)  we  have  for  the  steady  state  weight 


W^t) 


*  *  *  * 
KUq^  -IrI1  +S0S1  -SrS1 


(9) 


l+K(|lir+|S1, 


Substituting  in  (3)  and  using  (1)  we  obtain  the  following  for  the  output  of 
the  array  (as  derived  in  the  appendix  B-l,  (P-3)); 
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I0(i+g|s1|2)-gi1(s0s1i'  )+g(ir!i1|2+i1srs* ) 
i+g(|i1l2+|s1l2) 

Spq+glijVgs^iQi*  )^g(sr|s1l2^s1iri1* ) 
i+g(|i1l2+|s1l2 


+  N, 


gN1(I0I*  -IrI*  +sQs*  -SrS* 


i+gdiJ' 


'+  s. 


where  we  defined 

(1+Ko2)/K  -  1/g  -  Pth 


(10) 


(11) 


Pcli  stands  for  the  system  threshold  power.  The  first  term  of  (11) 
represents  the  contribution  of  the  interference  signal  I  at  the  output, 
while  the  second  is  that  of  the  desired  signal  S  and  the  last  terms  represent 
the  contribution  of  the  white  noise  at  the  input  to  the  output  of  the  array. 


The  Reference  Signal  Loop  Output 

The  output  of  the  reference  signal  loop  can  be  obtained  by  convoluting 
V(t)  with  the  impulse  response  of  the  loop  filter.  Therefore,  the  contibution 
of  the  desired  signal  at  this  output  can  be  written  as 

S  (l+g|l jVgS.U  I*  )+gS  (|S  |2)+gS1(I  I*  ) 

Sr  -  Fs  ± - £ - - - --S  1  (12) 

l+gdi^+lsj2) 

where  Fg  is  a  linear  operator  representing  the  effect  of  the  reference  loop 
filter  on  the  desired  signal  spectral  content.  When,  for  example,  the  local 
code  is  in  synchronism  with  the  incoming  desired  signal  code,  so  that  the 
desired  signal  has  a  very  narrow  band,  then  Fg  can  be  approximately 
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represented  as  a  constant  complex  number 


j0 

F  -  ye  Sr  (13) 

s 

where  y  and  0gr  are  the  gain  and  phase  shift  of  the  loop  filter  at  the 
desired  signal  center  frequency.  Combining  terms  we  can  write  (12)  as 

[l+g(!l1|2+|S1|2(l-Fs))]Sr  -  Fs[S0(l+g|l1|2)-gS1(I0I*1)+gS1(IrI1*)]  (14) 

Similarly,  for  the  interference  signal  we  get 

[l+g(|S1|2+|l1|2(l-FI))]Ir  -  F  [I0(l+g  S1  2)-gl1(S0S]*  )+gI1(SrS*  )]  (15) 

Equations  (14)  and  (15)  can  be  written  in  matrix  notation  as 


where , 


i+g(|s1|2+|i1!2(i-FI))  sf.V*-* 

I 

IftFj^S*  l+g(!l1l2+!si,2(l-Fs))  | 


(17) 


(18) 
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(19) 


isor-io(1+*IsiI2)-sii(sosi*)’ 

SNOR  ‘  V1<*II1|J)-  ‘Wl  >’  <20> 

F  0 

s 

B  -  ,  (21) 

-  0  Fj 

and  from  (P-6) 

A  -  [l+gd^^+lsj2)]  [l+g(|si|2(l-Fs)+|l1|2(l-FI)]  (22) 


Notice  that  when  B  *  0  the  reference  signal  loop  output  becomes  zero.  That 
is  to  say,  the  loop  filter  is  such  that  no  output  is  produced.  This  is 
obviously  the  case  when  the  filter  bandwidth  is  small  in  comparison  to 
both  the  desired  signal  and  interference  signal  bandwidths. 
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wx(t) 


-g  ( 


*  * 

IoIi  +sosi 

l+g(|l,|2+|S,]2) 


1S*  h  1  B  SoR 


NOR 


)  (24) 


where  S.,„  and  I  are  defined  in  (19)  and  (20)  respectively.  When  the 
NOR  NOR 

reference  loop  is  ineffective  (B  ■  0)  then 


'WT> 


-g  <V>o»1l>) 

W»(|I1|2+|S,|2 


(25) 


which  is  proportional  to  the  input's  cross  correlation.  However,  for  the 
general  case,  the  steady  state  of  the  complex  weight  becomes  (see  the 
derivation  of  P-8  in  Appendix  B-4). 


w 

where  we  notice  the  effect  of  the  reference  loop  filter  responses  (Fj  and  Fg) 
to  the  interference  and  desired  signal  respectively.  Particularly,  we  notice 
the  case  when  this  filter  has  perfect  response  to  the  desired  signal  (Fg  *1) ; 
then  the  weight  depends  on  the  cross  correlation  and  power  of  the  inter¬ 
ference  signal  only  and  hence  the  array  processor  will  produce  a  better 
null  in  the  direction  of  this  signal. 


(t) 


-g[IQIl  (1-F^SqSi  (1-Fa)] 
l+g[|S1|2(l-Fs)+|l1|2(l-FI/] 


(26) 


Desired  and  Interference  Signals  at  the  Array  Output 

To  obtain  the  signal  terms  at  the  output  of  the  array  we  can  use  (10) , 
in  which  we  substitute  the  expression  for  I  and  Sf  as  they  were  obtained  in 
(17).  However,  we  can  directly  use  (3)  to  get 


vs(t)  -  s0  +  Slwl(t) 


Vj(t)  «  I0  +  S^U), 


(27) 

(28) 


B-ll 


for  che  desired  signal  and  the  interference  signal,  respectively,  at  the 
output  of  the  array.  W  (t)  is  given  by  (25).  In  appendix  B- 5  we  derived 
the  final  expression  for  these  signals. 


ygd-FiM^liil^idpi;  )] 

l+gdsJ^l-F^  +  llJ^l-F^) 

Io+g(1-Fs)[Io!s1!2-Ii(sosi* )] 

l+g(|S1|2(l-Fs)+|l1|2{l-FI)) 


(29) 


(30) 


In  fact,  (29)  and  (30)  represent  the  complex  envelopes  of  the  desired  signal 
and  the  interference  signal,  respectively,  at  the  output  of  the  array  as  a 
function  of  the  complex  envelopes  of  these  signals  at  the  inputs  to  the  two 
elements  as  well  as  the  reference  loop  filter's  responses  to  these  signals, 

F  and  F  ,  respectively.  In  particular,  we  notice  that  when  these  responses 

S  x 

are  small  (Fg  1  0,  F^  -  0),  the  corresponding  array  outDut  signals  are 


VSN0R(C) 


i+g(!s1l2+|i1!2 


(31) 


V  (t) _ -225 _ 

IN0R()  i+g(|si|^|l1|2 


(32) 


We  will  further  assume,  without  loss  of  generality,  that  the  two 
array  elements  are  sufficiently  close  so  that  |l^|  ■  |lg|  *  I  and 
| Si]  ”  | Sq]  ■  S.  Together  with  the  definition  made  in  (2)  we  can  write 
(29)  and  (30)  as 
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(33) 


DESIRED  S IGNAL-TO- INTERFERED E  RATIO  AT  ARRAY  OUTPUT 

Since  the  different  terns  at  the  array  output,  (namely,  Vg(t), 

Vj(t)  and  V^t))  are  uncorrelated,  the  power  output  is  composed  of  the  sum 
of  the  powers  included  in  these  terms.  To  find  the  power  ratios  of  these 
different  terms  we  will  distinguish  between  different  processor  modes  of 
operation  by  using  the  different  ratios  of  the  bandwidths  of  desired 
signal  and  interference  signal  (at  the  input  of  reference  loop  filter) 
to  the  bandwidth  of  this  filter,  respectively.  Ideally,  the  reference  loop 


p 


filter  bandwidth  should  be  much  greater  than  the  desired  information 
bandwidth  and  much  smaller  chan  the  spread  signal  bandwidth. 


Acquisition  Mode 

Due  to  the  fact  that  the  local  spread  spectrum  code  is  not  in 
synchronism  with  the  incoming  desired  signal  code,  the  bandwidth  of  both 
the  desired  signal  and  the  interference  signal  are  larger  than  the  loop 
filter  bandwidth.  Therefore,  we  can  reasonably  assume  that  both  and  Fg 
are  approximately  zero.  Using  this  in  (33)  and  (34)  we  get 


| S 1 2out 

I  T  '  2 

I  I j  out 


•1(es  -0s 

H-alll2<l-e  10  10)12 


(36) 


,  J(6I  -9s  >, 

:n2|i+g!s|2(i-e  10  10)|2 


Computing  the  absolute  value  in  (36)  we  have 


si  2 out  *  Is!2 


l|2out 


i+2g|iru+g|in[i-cos(ec  -e  )] 

_ lo  no 

l[2  l+2g|s|2(l+g|s|2)[l-cos(9T  -0  )] 

lO  S10 


(37) 


Notice  chat  if  the  directions  of  the  two  signals  are  close  (cos(6  -6  )  *1), 

s10  L10 

i  i  2 

of  if  the  desired  signal  and  the  interference  signal  powers  (namely,  |S|  and 

t  i  2 

j 1 |  ,  respectively)  are  small  compared  to  the  threshold  power  Pth  (Pth  ■  1/g) , 
then 


out 


(38) 


'  ill  out  |I|- 

and  the  array  has  no  effect  on  the  output  signal-to-interf erence  ratio.  On  the 
other  hand,  if  the  relative  directions  and  the  powers  of  the  two  signals  are  such 
that  the  second  terms  of  both  numerator  and  denominator  of  (37)  are  much 
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larger  chan  one,  Chen 


1  s] 2quc  x  Is!2  1 i] 2(i+gl i] 2) 

!i!2ouc  !  i  ! 2  i  s  j 2  i-  i+g ;  s ! 2 ) 

5  111!  (39) 

!  s| 2 

This  is,  in  face,  che  pover  inversion  foraula. 

Thus,  Che  array  will  gradually  reduce  Che  power  of  Che  larger  signal 

•k 

which  is  necessarily  being  Che  inCerference  signal  leading  Co  acquisicion 

of  che  desired  signal  code  and  synchronizacion .  To  conclude,  we  depicc  in 

Figure  3  che  ouepue  signal-co-incerf erence  racio  as  a  funecion  of  che 

signal-co-Chreshold  power  racio  as  well  as  of  che  signal's  relacive  direction 

9*9  -6  .  This  is  done  bv  considering  the  different  signals-co-chreshold 

10  s10 

possibilities  from  (37), 


2.  Code  Tracking  Mode 

Due  to  synchronization,  the  bandwidth  of  the  desired  signal  at  the 
input  of  the  reference  loop  filter  is  small  (and  equal  to  the  information 
bandwidth)  and  can,  therefore,  be  assumed  to  be  much  smaller  chan  the  filter 


jes 

bandwidth.  Thus,  Fg  can  be  approximated  by  a  complex  number,  Fg  *  ye  r 
(see  (13)).  However,  che  interference  signal  bandwidth  at  the  input  of  the 
filter  is  assumed  to  be  much  larger  than  the  filter's  bandwidth  as  a 
consequence  of  the  spreading  caused  by  che  local  code.  Thus,  F^  -  0. 

Hence,  from  (33)  and  (34)  we  have 


* 

This  is  so  since,  if  the  desired  signal  were  the  larger,  then  the  svstem 
would  have  been  in  svnehronization  and  not  in  acquisition  mode,  as  it  is 
assumed. 


'  S 1 2out 


j(6s  6I  )  2 

.  sio  no  J 

Ugli;  (l-e _ It 


7  7  j®  -)  j(®T  )  I  2 

[I;  out  j  1 1  1  l+g(l-ye  Sr)|s',  (l-e  10  S10  )! 


Using  (p-11)  of  the  appendix,  this  becomes 


; S ) "out 

i  Ti  2 

, I ;  OUt 


,  (i+2g|i;2(i+gii;2)(i-cos(ec  -a.  ))] 

S|Z  S10  xi0 


I; 


l+2g|s|2f1+2g2js|4f2 


(40) 


where 


f.  *  1-ycos  6  -cos(0T  -8  )  +  ycos(9  +9T  -9  )  (41) 

1  Sr  ri0  s10  Sr  ^0  s10 

f.  m  (l+y2-2ycos9  )(l-cos(6  -9  ))  (42) 

Sr  Ho  S10 

First  we  will  consider  the  case  of  an  ideal  loop  filter:  that  is, 
when  y  -  1  and  Sg  =  0.  With  this,  f^  ■  f2  *  0  and  (40)  becomes 


2  2 

LS1.  ,  Ml.  [l+2g!  1 1  2(l+g|  lj  2)  (1-cos (9  -9  ))]  (43) 

( 1 1  ou t  |ir  io  ho 

Again,  if  the  angular  separation  of  the  two  signals  is  small,  or  if  the 
interference  signal  power  is  small  in  comparison  to  the  threshold  power,  then 

2  1 2  2  f  2 

| si  out/jlj  out  «  js(  / [ I j  and  the  array  has  no  effect  on  the  output  signal- 
to- interference  ratio.  Notice  that  this  result  is  true  regardless  of  the 
desired  signal  power.  If,  on  the  other  hand,  the  angular  separation  between 
signals  and  the  interference  signal  power  are  such  that  the  second  term  of 
(43)  is  much  larger  than  one,  then 


j  si  2 out  _  is;2  2g 2 [ l| 4(1-cos(9 

{  i  ! 2  ou  t  h!2 

For  this  to  hold  we  must  have 

1112  »  l/2g  -  Pth/2 


(44) 
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2  *2(2 

This  is  so  because  if  2gjl|  <<  1,  then  2g|l|  (l-t-gjl|  )  <<  1  in  contradiction 
to  our  assumption.  Notice  that  the  imorovement  in  the  desired  signal-to- 
interference  ratio  depends  on  the  square  of  the  interference-to-threshold 
power  ratio  ((g|l|  )  ■  (|l|  /Pth)  )as  well  as  on  the  angular  separation 
between  the  desired  signal  and  interference.  In  contrast  to  the  case  without 
reference  signal  (without  a  loop  for  extracting  the  desired  signal  or  when 
this  loop  is  ineffective  as  in  the  acquisition  mode  discussed  previously), 
this  result  is  independent  of  the  desired  signal  power.  Figure  4  depicts  the 
output  signal-to-interf erence  ratio,  for  the  ideal  filter  case,  as  a  function 
of  the  signal-to-threshold ,  the  interference-to-threshold  and  9  as  it  is 
described  in  (43). 

To  demonstrate  the  effect  of  a  nonideal  filter  in  the  reference  signal 

loop  we  reconsider  (40).  We  will  assume  that  e  <<  9  -9  so  that 

Sr  Ho  s10 

f.  -  (1-ycosS  )  (1-cos(9t  -5  ))  (46) 

1  sr  X10  S10 

With  this  (49)  becomes 

!.Sl2°u£  , 

1 1 1  Hut 

,  l+2g!  il  Hl+g|  lj  2)  (l-Cos(6  -9  T  )) 

I  s  1  _ Ho  Ho _ 

! X ! 2  l+[2g| s! 2(l-ycose  )+2g2:s|4(l-h-2-2ycosee  ) J (1-cos (r,  -7  ))  (4-) 

Sr  r  *10  S10 

In  order  to  minimize  the  effect  of  the  imperfect  estimation  of  the  reference 
signal  and  to  obtain  a  result  approximating  (43),  it  is  sufficient  to  recuire 
that  the  ednominator  factor  in  brackets  in  (47)  satisfy 
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(48) 


2g ! S | “ (l-ycosS  )+2g| s| ^(l+Y^-2ycos9  )  <<  1  (48) 

r  r 

Now,  let  us  further  assume  that  9g  is  very  small:  This  condition  might  be  met 


r 

by  using  a  reference  Iood  filter  wide  enough  for  the  desired  signal  information 
bandwidth  (but  narrow  enough  to  enable  the  assumption  of  F^  -  0) ,  or  by 
adding  (in  the  reference  signal  loop)  compensation  weight  that  will  reduce 


9s  [3].  With  this  condition  on  8  (48)  can  be  approximated  by 

r 

2g| S| “(l-ycos9  )+2g“ | S j ^ (1-ycose  )“  <  1/2 
r  r 

which  implies  that 


2  g : S |  (1-ycos  g  )  -  2 


or,  equivalently 

-  1.207  -  (1-Ycos  9  )  -  0.207  (49) 

;s'2  sr  'S'2 

Therefore,  for  larger  desired  signal-to-threshold  power  ratio,  vcos  9  must 

r 

be  closer  to  unity.  This  is  easy  to  establish  since  with  high  SNR  a  better 

estimator  for  y  and  9  can  be  obtained .  (For  example,  with  the  second 

sr 

compensation  scheme  of  [3]  it  was  shown  that  y  -*1,  9  "0,  and  obviouslv 

sr 

high  SNR  will  be  helpful.) 


The  inequality  condition  of  (48)  is  to  be  understood  as  searching  for  the 
assymptotic  behavior  of  the  denominator  of  (47).  Therefore,  for  3db  point, 
for  example,  it  is  enough  to  require  s  1/2  instead  of  <<  1  condition. 


[3]  V,  3ar-Ness,  "On  the  Problem  of  the  Reference  Loop  Phase  Shift  in  an 
N-Element  Adaptive  Array,"  l'P-VFRC-11-80 ,  July  1980. 
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Equation  (49)  represents  the  design  requirement  on  the  loop  filter 
response.  If  this  can  not  be  met  with  a  fixed  filter  design  then  a 
compensation  scheme  is  needed.  Notice,  however,  that  this  condition  is 
only  sufficient  (but  not  necessary)  to  obtain  the  desired  signal-to- 
interference  improvement  of  (43) .  Violating  this  condition  may  still  lead 
to  an  improvement,  although  of  smaller  value. 

If,  however,  (48)  can  r.ot  be  satisfied,  for  example,  because  the 
desired  signal  power  is  very  large,  then  we  might  have  instead  the  condition 
that 

g| S| 2(l-ycos0  )(l-cos(8_  -0  ))  >>  1  (50) 

sr  x10  s10 

Then  (as  we  show  in  (P-12)  of  the  appendix  B7),  the  denominator  of  (47)  can 
be  approximated  by 

2g2| S|4(l+y2-2ycos0  )(l-cos(0T  -8  )). 

sr  ■L10  s10 

Substituting  in  (47)  we  get 

LgJ2out  .Mi  i+2g|ii2(i+g!il2)  (51) 

1 1 1  2out  !  1 1  2  2g2|s|Zl(l+y2-sycos8  ) 

r 

2 

If  also  we  have  jlj  >>  1/g  ■  Pth  then,  after  some  easy  manipulation,  (51) 
becomes 

|sl2out  a  1 1!  ?  _ 1 _  (52) 

|l|2out  |S!2  (l-ycos0  )2+y2(l-cos@  ) 

sr  sr 

Notice  that  the  term  in  the  denominator  of  the  second  fraction  is  very  small 

(ycos0  -  1,  although  it  can  not  be  equal  to  one,  since  otherwise  (50)  can  not 
8r 

be  satisfied).  Comparing  with  (39)  we  observe  that  in  this  case  of  code 
tracking  mode  the  power  inversion  formula,  if  it  is  enforced  by  the  signal's 
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2  .  2 

input  conditions  (both  J I |  and  |s|  are  very  large),  is  still  favoring  the 
desired  signal. 


DESIRED  SIGNAL-TO-NOISE  RATIO  AT  ARRAY  OUTPUT 


First,  from  (35)  we  can  write  for  the  output  noise  power 


I  Vc)! 


(i  |  g2li0ii*  u-rp+SpS*  (i-fs)|2 
| l+g(| sl 2(1-Fs)+| Ix| 2(1-Fi) I 


where  we  used  the  fact  that  the  noise  processes  at  the  two  inputs  are 
uncorrelated.  Using  the  analysis  that  led  to  (P-13)  of  the  appendix,  we  get 


ai+a2 


|i+g[|sr(i-F#)+|iru-Fj)i 


(53) 


where 


[1+g 


Al+A2  -  l+2g[|srRel(l-Fs)+|l1|ZRel(l-FI)]. 
S|4!l-Fsl2+|l1!4|l-FI!2+|l1!2|si|2Rei((l+eJ(ai10  •H>>  )  (1-Fg)  (1-F*  ))] 


|  S | 2Rel[l-F8)+] I1| 2Rel (1-Fj) 


(54) 


Again  we  will  distinguish  between  the  different  modes  of  the  system's  operation; 
that  is,  the  acquisition  and  the  tracking  modes. 


1.  Acquisition  Mode 

Arguing  as  before  we  will  assume  in  this  case  that  both  F  and  F  are 

X  s 

approximately  zero,  so  that  we  can  write,  using  (33)  and  (53),  for  the  output 
desired  signal-to-noise  ratio, 
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N  out 


,  i+2g|i|2(i+g! i|2)(i-cos(ee  -e  )) 

|sr  _ 10  no _ 

o  2  l+2g(|s|2+!l|2)[l+g(|s|4+|l|4+|l|2|s|2(l+cos(6  -9  ) ) / I S | 2+ | I | 2 ) ] 

N  i10  sio 

2  i+2g|i|2(i+g|i|2)(i-cos(e=  -eT  ). 

>  S  _ s10  10  (55) 

aN2  l+2g( | S j  2+| l| 2) [ l+g( | S | 2+! 1 1 2) ] 


Therefore,  when  Che  directions  of  the  signals  are  so  close  (cos0_  -9  -  1) 

IlO  s10 

and  the  signal  power  is  small  compared  to  1/g,  or  both  signals'  powers  are 

2  2 

small,  then  the  output  SNR  is  bounded  below  by  |s(  /o^  which  is  the  input 
SNR.  When  the  angular  separation  between  the  two  signals  and  their  powers 
are  such  that  the  second  terms  of  both  the  numerator  and  the  denominator  of 


the  right  hand  side  of  (55)  are  much  larger  than  unity,  then 


(lsi2+|i|2)2 


dill 

2 

for  1 1 

°N 

=  Is!2  .  Iri4 

for  |  S 

That  is,  when  the  desired  signal  is  larger  than  the  interference  signal  at  the 
input  then  the  array  will  not  only  null  the  first  to  a  level  below  that  of  the 
second  (see  (39)),  but  it  will  also  degrade  the  signal-to-noise  ratio.  In  all 
other  cases  it  is  reasonable  to  expect  approximately  the  same  value  of  SNR  at  the 
output  of  the  array  as  at  the  input,  or  even  a  slight  improvement  in  SNR. 
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2.  Code  Tracking  Mode 


Again,  as  before,  we  assume  Fg  -  ye  r  and  F^  -  0  so  that  we  can 


write  for  this  case,  using  (33)  and  (53), 


! Sl 2out 
N  out 


,  i+2g| i| 2(i+g! il 2) (i-cos(e  -e.  )) 

1st _ sio  Ho _ 

crN2  l+2g(!l|2+!s|2  h1)+2g2(|l|4+|s!4h2+|s]2|l|2h3) 


where 


h^  «  1-ycos6s 


h  ”  (1+y  -2ycos6  ) 
i  sr 

h-  -  l-ycos9  +cos(0T  -0  )-ycos(0_  -0  +0  ) 


r10  sl0 


IlO  s10  sr 


For  the  case  of  an  ideal  reference  loop  filter,  0  ■  0,  y  *  1,  so 

sr 


that  hj  *  h^  ■  h^  *  0,  and  (57)  becomes 


,  ,  l+2g| I | 2(l+g| l| 2) (1-cos (0  -0T  )) 

l2out  _  _ S10  *10 


N  out 


i+2g|ira+g|i| 


Therefore,  except  when  the  directions  of  arrival  are  such  that  cos(0  -0T  )*1, 

s10  x10 

the  output  SNR  is  approximately  equal  to  the  input  SNR  regardless  of  the 
desired  signal  power. 

To  see  the  effect  of  a  non-ideal  filter  on  the  output  SNR  we  first  assume 
(48)  is  satisfied  and  write  (57)  as 


,  ,  l+2gil|2(l+g!l|2)a-cos(0T  -0  )) 

S  2ouC  .  isif _ _ _ 110  Sl0 

N  out  _  2  -  ~  ~ 

n  i+2g!i|2[i+g(!i'2+|s!2h3)] 

,  l+2g|  l|  ^(l+g(!  i|  2)  (l-COs(9  -0  )) 

.  1  si 2  _ X1Q  S10 

o  2  l+2g  I  2(l+g(  I  2+  S  2(1-  ) cos (  .  -  s  >) 

N  IlO  s10 
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<3 

where  bv  (48)  we  have  used  g| S 1 ‘ (l-ycos0  )  >>  1  in  h,  of  (38).  Coroaring 

sr  J 

with  (59)  we  notice  the  effect  of  y  i  1  on  the  output  SNR. 


Effect  of  Non-Zero  Interference  Residue  at  the  Output  of  the  Reference  Loop 

In  the  previous  discussion  we  assumed  that  in  the  tracking  mode  the 
reference  loop  filter  rejects  the  interference  signal  totally:  that  is,  we 
assumed  the  F^  -  0.  In  practice  this  will  not  be  the  case,  since  the  filter 
bandwidth  can  not  be  negligibly  small  (otherwise  it  will  reject  the  desired 
signal  as  well).  Thus,  •  ven  though  the  filter  bandwidth  is  smaller  than 
the  spread  interference  bandwidth  and  it  will  reject  a  large  part  of  the 
interference  power,  nevertheless,  there  will  still  remain  some  residue  at 
the  output  of  the  reference  loop  filter.  In  this  section  we  estimate  the 
effect  of  this  residue. 

We  will  assume  that  the  filter  is  ideal  in  estimating  the  desired 
signal,  that  is  Fg  »  1,  ,  however,  is  not  zero.  With  these  assumptions 

used  in  (33)  and  (34)  we  get 


I  S | 2out 
!  I i 2out 


,  j  (9.  -6  )  , 

l+g(l-FI|l|2(l-e  110  *10  )| 


Using  (P-11)  of  the  appendix  this  becomes 


[l+2gl l| 2e1+2g2| I 1 4e2 ) 


where,  analogous  to  (41)  and  (42),  we  define 

e,  “  1-8cos8t  -cos(8t  -6  )+Scos(8.  +eT  -8  ) 

1  xr  IlO  no  Ir  Iio  s10 

e,  -  (l+S2-2Bcos8T  ) (l-cos(6.  -8  )) 

2  lr  l10  s10 


(61) 


(62) 

(63) 
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r  '  « 


j9I 

Here  we  used  the  relation  1^(0  =  F^V^(t)  =  Se  rV^.(t)  where,  as  before, 

I„( t :  and  V^(t)  are  the  interference  signal  at  the  output  and  the  input, 

respectively,  of  the  reference  loop.  Note  that  V^(t)  is  also  the  array 

output.  Since  the  interference  signal  is  assumed  to  be  wideband  (it  is 

widened  by  multiplication  with  the  system  code),  6  and  9  would  obviouslv 

xr 

be  frequency-dependent.  As  a  first  approximation  we  may  assume  that  the 
reference  loop  filter  has  an  ideal  rectangular  shape  with  zero  passband  phase 
shift.  That  is,  we  assume  0T  -  0  and  8*1  within  the  passband  and 
3*0,  and  9^  is  undefined  outside  the  passband. 

Define 


S  -  (l/BI)/9(w)dw 


(64) 


where  is  the  interference  signal  bandwidth  and  B^  is  the  reference  loop 
filter  bandwidth.  Then  (61)  together  with  (62)  and  (63)  becomes 


2 

Ili°.ut  .  Is \  [i+2gll!2(l-8)(l+g[ll2(l-B))(l-cos(eT  -0  ))] 

| I | out  |I|“  *10  s10 

Comparing  with  (43)  which  was  obtained  with  an  assumption  of  an  ideal 
filter  (F  *  1),  we  notice  that  the  effect  of  the  interference  residue  at 
the  output  of  the  reference  loop  is  equivalent  to  that  of  reducing  this 
signal  power  by  1  -8. 
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Figure  2  The  LMS  Interference  Canceller  with  Coded  Reference 
Signal  Loop  for  Spread  Spectrum  Communication 


-->6 
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APPENDIX 


B-l  The  arrav  output 


v(t)  =  XQ(t)+W1(t)X1(t)  (P-1) 

Using  the  definition  of  X^Ct),  i  =  0,1  together  with  W^lt)  from  (6)  and 
corresponding  terms  of  (7)  and  (8),  we  get 

*  *  *  * 

R(I.I.  -II.,  +S.S.  -S  S  ) 

V(t>  =  I0+So+No - - J-± -  (I,+S,+N,)  (P-2) 


l+Kd^l^+lS^+o^  ) 


1  "1  ‘ V 


Let 


(1+Kc^  / K  =  1/g  -  Pth 


wherfe  Pen, stands  for  Che  system  threshold  power.  Then  (P-2)  becomes 

8(I0Ii*"Irll  +S0Si*-SrSi*) 

V ( t )  -  (In+S0+N0) - ^ - — - - - -  (I ,+S .+N  ) 

0  0  0  i+g(!i  ;2+>s  ,2)  1  1  1 


(Vs0)ti>g(!i1l^!sir)i-g(Il+s1)[i0i1  -iri1^s0s1  -srsL 

1+8(1 IL| 2+i sL! 2) 


+V 


Sha0ll-lrh+50Sl  -SrSl  } 


'  t  !  2-i-!  s  I  2 
l1 


l+gd^r+is 


n-  )o 


Thus 


V(t) 


i0(i+gl sj 2)-gi1(s0s*  )+g(Irl ^ 
i+g(|s1!2+|i1!2) 


B-3  Using  (18)  we  have 


*  ,  *  -1 
[S1  !  l1  ]  A 


s*  [i+g(|s1|2+U1!2(i-FI)] 


*  i  1 2 


*  I  1 2 

Fs|Sl>  + 


I*  (l+g(jl1l2+iS1|2(l-Fs) 


i+g(!s  |2+I,I  I2)  *  .  * 

=  - —  is  •  \ 

5-4  From  (24)  and  using  (21)  we  write 


W,  (t)  *=  -g(- 


(  V1V1* 

l+s(l  IjJ  2+i  Sj 


FsSl  smor+fiii  xnor. 


Using  the  definitions  of  S^0R  and  1^0R  in  (19)  and  (20),  resDectively ,  and 
that  of  4  in  (22)  we  can  write 


WL(t)  =  =*-  u0i*+sQs*  )[l+g(!si!2(l-Fs)+;ii;2(l-FI))] 


Vi  [S0(i-8ji1!2)-gs1(x0x‘ )] 


FIX1  [I0(1+g'Sll  "8ll(S0Sl  )]} 


■  ~f~  ' IQIi*  [l+g(iS1i2(l-Fs)+!l1!2(l-FI))+gFs!si]2-FI(l+giS,|2;] 


S0S*  [l+g(:S112(l-Fs)+!li;2(l-FI))+gFI|l1!2-Fs(l+g!l1|2)]: 
which  can  be  brought,  after  some  simple  manipulation  and  (22),  to  the  form. 


-glio^O 

l+e[ '  S  ;  2 <"  1— F  )-!-',  I,  2 ( 1— F _ )  ] 
1  si  I 


B-5  Substituting  (26)  in  (27)  we  have 


*1 


„  m  -  .  ^O1!  (1-FI)4S0S1  “"V1 

8  '  0  l+g(|S1!2(l-Fs)+|l1|2(l-FI)) 

[l+g(|S1!2(l-Fs)+!l1|2(l-FI)]S0 
l+g(lsi!2(l-Fs)+|l1|2(l-FI) 

^Wl  (l-FI)-giS1|2(l-Fs)S0 
l+g(|S1|2(l-Fg)+|l1|2(l-F1) 

_  S0+8(l-FI)[S0ll1l2-S1(I0Ii*  >] 

l+g(|S1|2(l-Fs)+!l1|2(l-FI) 
and,  similarly,  we  obtain  by  using  (26)  in  (28) 

. 

v  (t)  -  - = - = - 

i+g(|s1r(i-Fg)+|s1ra-F1)) 

B-6  Let  g(a,F ,0)  -  l+a(l-F) (l-ej6 ) 
where  F  is  complex  and  6  are  reals.  Then, 

j  g(a,F,0)|  -  [1+  (l-F)(l-eJ0)][l+  (1-F*  (l-e“j0)] 

-  1+2  (l-Rel(F)-cose+Rel(Fej9)] 

+2  2 (1+| F| 2-2Rel (F) ) (l-cos6 )  (P-11) 

I  I  ° 

B-7 ^  Condition  (50)  obviously  implies  that  g|Sj  (1-  cos9  )  >>  1. 

SY* 


(P-9) 


(P-10) 
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(1+v  -2ycos9  ) 
sr 


l-ycos9 


g \ S ' " (1-ycose 


1+y  -2ycos9 


Sr  (l-ycos9 s  )2 


>  g  S  “(l-ycos9  )  >>  1 


Therefore,  the  denominator  of  (47)  can  be  approximated  by 


( l+y*-2y  cos6  ) 

i  2  2  sr 

l+2g  S!  (l-vcos6  )[l+g!S!  - ](l-cos(9  -9s]_q^ 

Sr  1-y cos9 


-  2g2'  S'.  4(l-^v2-2vcose  )  (l-cos(9  -9  ))  (P-12) 

sr  1 10  s10 


3-8  To  derive  the  noise  power  at  the  output  of  the  array  we  notice,  using 
(35) ,  that 


8‘Vl  (1-FI)+S0S1  <1”Fs>!‘ 


'v  (t)  '  -  ?  ‘  (1+  - V-i - i - - - ^-r- 

!i+g(|s1i2(i-rs)+|i1|2(i-F1))|2 

where  we  used  the  fact  that  N.  and  N  are  uncorrelated. 

0  1 


Aj_  i  jl+g(iS1|2(l-Fs)+|ll2(l-FI))|2 

-  l+g2[|S1|4|l-Fs!2^!l1!4!l-FIl2+2!l1!2,S1l2Rel(l-Fs)(l-FI) 


+  2g[ | S1( 2Rel(l-Fs)+| ll [ 2R«1(1-F1) I 
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ABSTRACT 

The  performances  of  two  Rinds  of  interference  cancelers  are 
compared,  namely,  an  interference  canceler  that  uti.  izes  a 
directional  constraint  (directionally  constrained  interference 
canceler)  and  an  interference  canceler  with  reference  signal  loop 
(LMS  interference  canceler).  For  simplicity,  this  is  done  for  a 
two-element  array.  In  our  comparison  we  use  the  array  output 
desired  signal-to-interf erence  power  ratio  (SIR)  and  desired  signal- 
to-noise  ratio.  These  power  ratios,  obtained  with  the  first  canceler 
when  the  angle  of  arrival  of  the  desired  signal  relative  to  the 
constraint  direction  is  sufficiently  small,  are  compared  with  the 
corresponding  power  ratios  obtained  with  the  second  canceler  when 
it  operates  in  code  tracking  mode.  Similarly,  power  ratios,  obtained 
with  the  first  canceler  when  the  angle  of  arrival  relative  to  the 
cdnstraint  direction  is  not  so  small,  are  compared  with  the  corres¬ 
ponding  power  ratios  obtained  with  the  second  cancclcr  when  it 
operates  in  code  acquisition  mode.  Finally,  in  considering  the 
"accurately  constrained"  case  of  the  first  canceler and  the  "ideally 
filtered"  case  of  the  second  canceler,  we  establish  the  condition 
under  which  one  or  the  other  of  these  cancelers  perfor-.s  better. 


* 


On  leave  of  absence  from  the  School  of  Engineering,  Tel  Aviv  University 


INTRODUCTION 


Adaptive  array  processors  are  used  more  and  more  as  a  method  for 
interference  cancellation  in  the  design  of  advanced  systems.  However, 
these  automatic  processors  might  also  cancel  the  desired  signal  unless 
it  is  somehow  protected  by  the  algorithm  used:  in  particular,  this  is 
true  when  the  desired  signal  power  is  much  larger  than  the  interference 
power.  One  way  of  obtaining  such  protection  is  bv  using  a  cross-coupled 
processor  that  implements,  in  a  genuine  boot-strapping  form,  the  "power- 
inversion"  propertv  of  the  interference  canceler.  This  form  of  "power 
separator"  was  reported  in  [1].  Customarily ,  some  a  priori  known 
properties  of  the  desired  signal  are  used — properties  which  sufficiently 
distinguish  it  from  other  signals  (interferences). 

An  array  for  spread-spectrum  communication  system  has  been  described 
[2,3]  using  a  reference  signal  loop  to  process  the  output  of  the  array 
and  extract  a  sufficiently  clean  estimate  of  the  desired  sign.,..  This 
is  then  used  as  a  reference  in  the  Widrow  LMS  algorithm  [4].  Such  an 
adaptive  array  processor,  as  it  was  proposed  by  Compton  [2],  is  depicted 
in  Figure  1.  In  Figure  2  a  different  approach  is  presented:  The  approach 
is  termed  "LMS  interference  canceler,"  and  its  performance  was  studied  in 
[3].  Both  approaches  use  the  special  spectral  property  (spread  spectrum) 
of  the  desired  signal  to  keep  the  processor  from  affecting  the  desired 
signal  regardless  of  its  dynamic  range.  As  a  result,  the  desired-signal- 


to- interference  ratio  (SIR)  is  substantially  improved  except  when  the 
relative  directions  of  arrival  of  these  two  signals  are  close.  The 
improv. ment  in  SIR  depends  on  the  signals'  relative  directions  of  arrival 
as  well  as  on  the  interf erence-to-threshold  power  ratio  (ITR) :  (Threshold 
power  approximately  equals  the  system  noise  power). 

The  improvement  decreases  when  the  estimation  of  desired  signal  by 
the  reference  signal  loop  becomes  less  and  less  accurate. 

An  alternative  use  of  a  priori  information  is  found  in  the  constrained 
algorithm  [5].  Here,  knowledge  of  the  desired  signals’  direction  of  arrival 
is  used  to  prevent  Che  processor  from  affecting  the  signal.  The  canceler 
of  this  approach  will  be  called  "directionally  constrained  interference 
canceler . " 

In  radar,  the  direction  of  arrival  of  the  desired  signal  is  assumed 
accurately  known  (the  direction  of  the  main  beam)  and,  therefore,  the 
second  alternative  can  be  used  satisfactorily.  In  point-to-point  communi¬ 
cation,  on  the  other  hand,  the  direction  of  arrival  of  the  desired  signal 
might  be  known  with  only  a  modest  level  of  accuracy.  Some  special  spectral 
property,  such  as  spread-spectrum  modulation,  might  or  might  not  be 
available.  Nevertheless,  the  questions  that  the  system  designer  might  be 
faced  with  are  (considering  the  output  SIR)  which  of  these  approaches  is 
preferred?  How  does  the  preference  depend  on  the  dynamic  range  of  the 
input  powers?  How  does  it  depend  on  uncertainty  in  the  desired  signal's 
direction  of  arrival  (used  in  the  second  approach).  How  does  it  depend  on 
the  amount  of  imperfection  of  estimating  the  reference  signal  (used  in  the 
first  approach)? 

The  performances  of  these  two  alternatives  are  to  be  compared  in  this 


article.  We  first  use  the  Applebaum  approach  for  a  directioanl  constraint 
interference  canceler  and  derive  the  output  desired-signal-to- interference 
power  ratio  (SIR)  and  desired  signal-to-noise  ratio  (SNR).  For  simplicity 
this  is  done  for  a  two-element  array.  In  particular,  we  notice  the  change 
in  performance  due  to  change  in  desired  signal  direction  relative  to 
constraint  direction.  We  then  compare  the  performance  of  the  directional 
constraint  canceler  (hereinafter  called  canceler  A)  with  that  of  an  inter¬ 
ference  canceler  that  utilizes  a  reference  signal  loop  (hereinafter  called 
canceler  B).  In  our  comparison  we  use  the  array  output  SIR  and  SNR.  These 
power  ratios,  obtained  with  canceler  A  when  the  angle  of  arrival  of  the 
desired  signal  relative  to  the  constraint  direction  is  sufficiently  small, 
are  comapred  with  the  corresponding  power  ratios  obtained  with  canceler  B 
when  it  operates  in  the  code  tracking  mode.  Similarly,  power  ratios 
obtained  with  canceler  A  when  the  angle  of  arrival  relative  to  the  constraint 
direction  is  not  so  small  are  comapred  with  the  corresponding  power  ratios 
obtained  with  canceler  B  when  it  operates  in  code  acquisition  mode.  Finally, 
in  considering  the  "accurately  constrained"  canceler  A  and  the  "ideally 
filtered"  canceler  B,  we  establish  the  conditions  under  which  one  or  the 
other  of  these  cancelers  performs  better. 

DIRECTIONALL  CONSTRAINED  INTERFERENCE  CANCELER 

The  array  arrangement  is  depicted  in  Figure  3.  In  complex  envelope 
notation  the  inputs  to  the  array  elements  are  denoted  by  X^,  X,,  ...  X^ 
where  n  is  the  number  of  elements  and, 

m 

x  (t)  -  s.(t)+s  (t)  +  y  i  .(t) 

j-1  J 


C  — 4 


Cl) 


I  (t) ,  j»l,2,...m  is  the  jth  interference  signal  at  the  array's  ith  element, 
Si(t)  is  the  desired  signal  at  the  ith  element  and  N^(t)  is  the  noise  term. 

In  vector  notation 

m 

X(t)  -  S(t)+N(t)+  l  I.(t)  (2) 

j-l° 

For  plane  wave  signals  we  can  write 

S(t)  -  s(t)S  (3) 

I.(t)  -  i  (t)Ij  (4) 

where  S  and  1^  are  the  complex  magnitude  vectors  (phase  and  amplitude)  of 
the  desired  signal  and  the  interfering  signals,  respectively.  These 
vectors  obviously  depend  on  the  direction  of  arrival  of  these  signals  with 
respect  to  the  array.  We  will  also  assume  that  | s(t ) | 2  =1  and  ] i ( t j [ 2  -1. 

A  processor  for  a  directaionally  constrained  interference  canceler  is 
one  that  derives  the  weight  vector  W,  such  that  in  the  steady  state  the 
following  statements  are  satisfied  simultaneously; 

|V(t)l2  -  I W^X ( t ) 1 2  «  min  (5) 

and 

E^W  *  n  (6) 

where  T  stands  for  transpose  and  the  overbar  for  the  mean  value.  That  is, 
according  to  (5)  the  processor  will  minimize  the  mean  square  value  of  the 
array  output.  3v  (6),  any  signal  s(t)  impinging  on  the  array  from  a 
direction  E  will  not  be  affected  by  the  array  (will  become  s(t)’n  at  the 
array  output;.  Therefore,  E  is  the  cLOizction  0&  coniticunt. 

c-r> 


For  the  implementation  of  this  processor  we  propose  to  use  the 
’raainbeam  constraint"  approach  of  Applebaum  and  Chapman  [5]  as  it  is  given 


in  Figure  4.  The  output  of  the  summer,  Y  ,  where 


YQ(t)  -  $  X(c)  (7) 

($  is  the  beam  steering  vector) ,  and  the  outputs  of  the  linear  transforma- 

t ion  Y ( t )  »{l(t),...,  Y  ,  ( t )  }  ,  where 
.  i  n-i 

Y(t)  -  AX ( t) ,  (8) 

are  the  main  and  auxiliary  inputs  of  an  unconstrained  interference  canceler. 
The  weight  vector  of  this  canceler  W  -  {W  }  is  controlled  so 

-y  yi  y 

n-1 

that  the  mean  square  value  of  the  arrav  output  |v(t)|  is  minimized.  That 
is,  in  the  steady  state  must  satisfy 

[l+KY*(t)YT(t)]W  «  -KY*(t)Y  (t)  (9) 

-y  -  ° 

where  K  is  the  gain  of  the  control,  loop,  I  is  an  (n-1)  identity  matrix 
and  the  (*)  stands  for  the  complex  conjugate. 

It  is  possible  to  show  [5]  that  if  we  choose  the  steering  vector 


d  -  E 


che  direction  of  constraint,  and  choose  the  transformation  A  such  that 

AE  «  0,  (11) 

then  Che  solution  of  (9)  together  with  (7)  and  (8)  is  the  solution  of  the 
constraint  problem  (5)  and  (6),  which  we  called  the  "directionally 
constrained  interference  canceler."  Notice,  from  Figure  4  and  (8),  that 


Hence,  by  (10)  and  (11), 


A£  =  0  =*  B1  -  0 

(where  1  =  (1,...,1)  is  Che  required  constraint  transformation. 


(13) 


The  Two-Element  Case 

To  facilitate  a  comparison  of  the  performance  of  the  directionally 
constrained  interference  canceler”  with  that  of  the  LMS  interference 
canceler,"  we  will  restrict  ourselves  to  the  case  of  a  two-element  array 
and  a  single  interference  signal.  That  is,  in  (2)  the  sum  contains  only 
one  element  and  the  vectors  are  two-dimensional.  Also,  without  loss  of 
generality  we  will  assume  that  the  direction  of  constraint  is  broadside 
to  the  array.  That  is,  from  (10)  $  -  E  *  1,  and  from  (12)  A  -  B.  To 
satisfy  (11)  we  must  use  A  -  [1,-1],  and  by  (7)  and  (8), 

Y  (t)  =  X, (t)+X-(t) 

O  1  * 

and 

Yx(t)  -  X1(t)-X2(t)  (17) 

where  X^t)  is  defined  by  (1)  with  m-1  (one  interference).  This  arrav  arranee- 
raent  is  given  in  Figure  5.  Obviously,  for  this  special  case  (9)  becomes 


[l+KY*(t)Y1(t)]Wy  -  -KY*U)YQ(t)  (18) 

it  ★ 

Using  (16),  (17)  and  (1)  in  the  derivation  of  Y1(t)YQ(t)  and  Y^OY^t) 
in  Appendix  A-l,  equation  (13)  can  be  written  as 

-g[(S1-S2)*(S1+S2)+(I1-I2)*(I1+l2)l 

w..  “  o 

i+g(is1-s2!  +|i1-i2r) 
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(19) 


! 


where  - 

,  a 

1/s - j - Pth  (20) 

stands  for  the  system  power  threshold-  and  S2,  1^  and  1^  are  the  com¬ 
ponents  of  vectors  S  and  I,  respectively,  as  they  are  defined  in  (3)  and  (4). 
The  Array  Output 

From  Figure  5  we  have 


V(t)  -  YQ(t)  +  Yx(t)Wy 

Using  (16)  and  (17)  with  (1)  we  show  that  (see  Appendix  A-2)  , 


(21) 


S  (t)[l+2g(|l  |2-I*I  )]  +  S7  ( t)  [  l+2g  ( 1 1..  |  2-I*I  )  ] 
V(t)  -  -i - 1 - - £ - - - i - L_L_ 

l+gdSj-S^  +  lix-i^r) 

I1(t)[l+2g(|s2|2-S*S2)]  +  I2(t)[l+2g(|si!2-S2S1)] 

I^g(l s1-s2l 2  +  i I1-X2I 2) 


H1(t)[l+2g(|S2|  + 


il. 


,2 


I1I2^  ^ 


N2(t)U+2g(iS1i2-S*S1  -,Ir‘--I2I1)] 


usd  s  s  2 ! ' 


1 1  -1  is 

ii  2  ; 


(22) 


Clearly,  the  first  term  is  the  desired  signal  at  the  output,  the  second  term 
is  the  interfering  signal,  and  the  last  term  is  the  noise.  We  will  further 
assume  that  the  two  array  elements  are  sufficiently  close  so  that 
|S1i  »  |s,|  -  1  s|  and  UJ  -  1 12 !  -  |l|.  Hence, 


S 

s 

I 

I 


1 

2 

1 

2 


is|e 
|3|e 
1 1 !  e 


*6s 

*6s 

*9i 


1 

1 

2 


i  23^ 
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where 


A 

ai  d 

< 

CD 

1 

a> 

-  6 

=  siniji 

Si 

s 

c  s 

9-1-9-0 

il  i2 

A 

-0. 

1 

(u.d 

=  -  sirnii . 

C  1 

(24) 


d  is  che  distance  between  the  array  elements,  i*)g  and  <u^  are  the  radian 
frequencies  of  the  desired  and  interference  signals,  respectively,  c  is 
the  speed  of  electromagnetic  waves  and  ij>s  and  are  the  directions  of 
arrival  of  these  signals  with  respect  to  the  array  broadside  direction, 
respectively.  With  this,  (22)  can  be  rewritten  as 


[S  (t)+S  (C)](l+2g!in  -  2g(I  I  S.(t)  +  I,I.S,(t)) 

V(t)  -  — ^ ^ - 5 - 1-2-1-, - 2_i_2 -  (25) 

i+sd  s1-s2i z  +  iiri2r) 

[I1(t>  +  X 2 ( t ) ] (l+2g| S| 2)  -  2g(S*S2I1(t)  +  S*S1I2(t)) 

i+gd  s1-s2! 2  +  ix1-x2;2) 


N1(t)[l+2g(|Si4  -  S"  S2+|lix-lJl2)]  +  N2(t)[l+2g(!sj2-S2  S^i'lj^-I^)] 

i+S(|s1-s2;2  +  i  2) 

Using  the  derivation  in  Appendix  A-3,  we  finally  write 


,2  * 


V(t) 


(26) 

j/2(9  +0 

2 s ( t )  ]  S  |  e  Si  s':',[(l+2g[lp)cos(6s/2)  -  2g|  I  j  cos(9  . -  9g/2 )  j 
l+2g[  ]  S]  2(l-cos9g)  +  !  I  ]  ^(1-cosS^  ) 
j  /2(9  +  6.  ) 

2i(t) 1 1 ; e  11  l2  [ (l+2g j S i 2)cos (5  /2)  -  2g j S i  2cos (9 g- 9./2 ) ] 
l+2g  [  1  S  j  2  ( l-cos9  s )  -t-  |  i]  2(l-cos0J  ] 

S1(t)[l+2g(.S  (1-e  V  i;“(l-e  1) )  ]+N’2  ( t)  [  l+2g  ( 1  S  ( 1-e  s)+!Ij“(l-e  1))] 

l+2g [ | S J 2 ( l-cos9  )  +  | I | 2 ( l-cos6  ) ] 

From  (26),  noting  that  the  three  terms  are  uncorrelated,  we  obtain 


09 


after  some  trigonometric  manipulation 


improvement  factor  becomes 


'  do 


fl+4g|r|,2sin2(0i/2)]2 
cos2(6 ,/2) 


(31) 


The  Desired  Signal-to-Noise  Power  Ratio  at  the  Arrav  Output  (SNR) 
From  (27)  and  (29)  we  get 


■  c,  2  „ 
i  S i  out 

Nout 


cos  (9  /2) [ l+2g  I  “sine . (tan(9 . /2)-can(9  /2 ))]' 
S  1  1  s 


l+4gh^+8g 


(32) 


h^  =  j S| 2(l-cos6s)+j l| 2(l-cos6^) 

h.  =  | Si 2| l| 2(1+cos(8 .-0  )-cos6  -cosB , )+| s| ^ (1-cosB  )+| I j ^ (l-cos9 . ) 
2  1  is  s  i  1  s  l 


PERFORMANCE  COMPARISON 

We  now  intend  to  compare  the  performance  of  the  directionally  constrained 
interference  canceler  with  the  performance  of  the  canceler  which  utilizes  a 
reference  signal  loop.  First  we  compare  the  desired-signal-to-interference 
power  ratio  of  these  two  cancelers  and  then  consider  their  desired-signal-to- 
noise  power  ratio 

Comparing  the  Desired-Signal-to-Interf erence  Power  Ratio 

In  doing  so  we  will  distinguish  between  two  cases,  where  the  desired 
signal  angle  9g  (relative  to  the  constraint  direction)  is  either  small  or 
not  small. 

For  the  second  case  we  may  notice  the  following  from  (30): 

(1)  If  the  relative  direction  between  the  desired  signal  and  inter- 


ference,  and  the  powers  of  the  two  signals,  are  such  that  the  second  terms 
in  brackets  of  both  the  numerator  and  denominator  of  (30)  are  much  larger 
than  one,  then  one  can  easily  show  that 


'  dl 


|  X|A  sin2(6i/2) 
1  Si 4  sin2 (8g/2) 


and  that 


s i ^  out  iii2 


out 


[  S I  2  sin2 (0  /2) 
s 


(33) 


(34) 


That  is,  besides  the  "power  inversion"  relation  we  have  further  improvement  in 
SIR  since  one  might  assume  0g  <  9^. 

(2)  If  the  signals'  directions  are  close  and/or  the  signal's  powers 
are  small,  such  that  the  second  terms  in  brackets  in  both  the  numerator  and 
the  denominator  of  (30)  are  much  smaller  than  one,  then 


cos~(9g/2) 
d2  cos2(9./2) 


(35) 


and  that 

■  Isli  -  (36) 

|l[^out  [ I 1  cos  (0^/2) 

That  is,  we  have  some  improvement  due  to  the  assumption  that  ?g <  5^. 

It  is  reasonable  to  compare  the  performance  of  the  directionally  constrained 
canceler  in  this  case  (0  is  not  very  small)  with  the  performance  of  the 
other  canceler  only  when  this  other  canceler  operates  in  the  acquisition 
mode.  That  is,  when  the  reference  loop  generates  no  reference  signal.  In 
so  doing  v®  first  recall  (see  [3])  that  under  the  condition  of  (1)  and  (2) 
above,  the  corresponding  output  desired-sienal-to-interference  power  ratios 


0 1  2 


r 


of  the  interference  canceler  with  reference  signal  are  given  by 


(1) 

(2) 


S  j  2out  _  1112 
1 1 2 out  ! S I  2 

sl 2out  _  | S I  2 
l|2out  |l| 2 


Thus,  this  last  canceler,  when  it  is  in  acquisition  mode,  is  inferior  to  the 
directionally  constrained  canceler  provided  the  direction  of  the  desired 
signal  is  closer  to  the  constraint  direction  than  the  direction  of  the 
interference  signal. 

To  consider  the  case  where  the  desired  signal  angle  of  arrival  9g 
(relative  to  the  constraint)  is  small,  we  first  write  (30)  as 


2g!l,|2sin9  can(9  /2)  ,  - 

[1 - ~~i- - 2 - J  2  cos  8  /2 

:i+4g!li  sin  (6/2)]*  (l+4g|ll  sin  (et/2)) 


cos  (S./2) 


l+4g ■ S | 2sin2 (e  /2)(l-tan(3./2)/tan(9/2))]^ 


In  oarticular,  if  4g | I | -sin^6^/2) is  much  larger  than  one,  then 


F  *  F 
d  ‘do 


cos2(?s/2)[l-tan(eg/2)/tan(3./2)]2 
[l+4g|S|2sin2(9s/2)(l-tanO ,/2)/tan(6s/2))]2 


(38) 


where  F.  is  the  improvement  factor  when  9  "0,  and  the  second  term  repre- 

00  s 

sent.*  the  change  due  to  a  non-zero  9g  (i.e.,  a  non-accurate  constraint). 

For  9  <9.,  and  if  9  is  sufficiently  small  so  that 
s  l  s 

sin2(9  /2)  +  — —  sine  <  l/4g | S I  2 ,  (39) 

s  2v~ 

then  the  improvement  factor  in  dB  is  given  by  (see  Appendix  A-4) 


oi  3 


T 


<rd>dB  "  201o«[' 


l+4g | I | “sin2 ( 0 . /2) 


co s 6  .  / 2 
i 


-10{ (es/2)2[l+8g!s!2]+2(9s/2)f l/tan(e./2)-4g!s!2tan(e./2)] }  (40) 


Notice  that  (F^)^  is  a  quadratic  function  of  eg/2  which  depends  on  the  loca¬ 
tion  of  the  interference  8^  and  the  desired-signal-to-threshold  power  ratio 

,  .  2 

(g | S | *)  ,  but  not  on  the  power  of  the  interference.  This  function  has  a 
maximum 


l+4g!lj2sin2(9./2)  [l/tan(e./2)-4g!s!2tan(9./2)] 

([Fd]dB)m  *  201o«[ - ~T7T~r~"~ —  1+10 - 1 - : - 1 - (41) 


cos(6i/2) 


l+8g  |  S 


:  2 


at 


(e  /2)  -  - 

s  m 


[l/tan(6i/2)-4g|srtan(e./2)] 
l+8g | S 1 2 


(42) 


That  is,  depending  on  whether  1 1 an e ^ / 2 |  is  greater  or  less  than  1/2  /g | S | 

(0  12)  has  the  same,  or  opposite  sign  as  that  of  9 . .  9  *  0  only  when  these 

s  m  rr  °  l  sm  * 


terms  are  equal. 


It  is  reasonable  to  compare  the  performance  of  the  directionally  constrained 
canceler  when  8g  is  small  with  the  performance  of  the  other  canceler  (with 
reference  signal  loop)  operating  in  the  tracking  mode.  Recalling  [3]  that 
the  corresponding  imorovement  factor  of  the  interference  canceler  with  reference 
signal  loop  is  given  by 


l+2g 1 1 1  2 (1+g 1 1 1 2) ( 1-cos ( 8-8 . ) ) 

F  -  - 5 -  - -  (43) 

l+2g|s|2f1+2g!s|4f2 

where 

f,  -  1-YCOS0  -cos(9  -9,)+ycos(9  +6-9  )  (44) 

i  sr  s  l  sr  i  s 

f-  *  (l+y2-2ycos0  ) (1- cos (9  -8  ))  (45) 

^  sr  1  s 
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is  Che  reference  loop  filter  complex  envelope  response  to  the  desired 


signal. 

suffers 

residue 

without 


That  is,  we  assumed  that  the  desired  signal's  complex  envelope 
phase  shift  5^  and  attenuation  y ,  while  the  interference  signal 
at  the  output  of  the  loop  is  small  and  can  be  ignored.  We  may  assume 
loss  of  generality  that  9@  =  0.  For  an  ideal  loop  filter  we  have 


Therefore , 


F 


R0 


is 


9 


sr 


<  < 


9.-6 
l  s 


l+2gjl|2(l+g|l|2)(l-cose.)  (46) 

9^  we  can  write  (43),  together  with  (44)  and  (45), 


F  =*F 
R  R0 


l+[2g | S | "(l-ycosesr)+2g  | S]  ■,(l+Y"-2ycos0sr) ] (l-cos9i) 


(47) 


The  second  factor  represents  the  change  in  the  improvement  factor  due  to  a  non¬ 
ideal  filter.  If  the  term  in  parenthesis  in  the  denominator  is  smaller  than 
one,  we  can  write  for  the  improvement  factor  (in  dB) , 


(VdB  °  101oRU+2g|l|2(l+g|l!2)(l-cosei)] 

-  20( [g | S | 2 (l-ycos9sr)+g2 | S I ^ (l+y2-2ycos9sr) ] (l-cose^) ] )  (48) 

Notice  that  (?_),_  is  a  decreasing  function  of  ycos9  (y<l)  and  has  its 
R  do  sr 

maximum  at  y«l ,  9gr“0  (ideal  filter).  The  degradation  in  performance  is 

2 

steeper  when  the  desired-signal-to-threshold  power  (*g|S|  )  is  larger. 

Finally,  we  compare  the  performance  of  the  two  cancelers  in  the  case 
when  the  first  is  accurately  constrained  and  the  second  has  an  ideal  filter. 
The  corresponding  improvement  factors  are  then  given  by  (31)  and  (46),  re¬ 
spectively.  Under  these  conditions,  the  first  has  better  desired-signal-to- 
interference  power  ratio  if  and  only  if 


C-15 


In  Appendix  A-5  we  show  that  for  the  validation  of  (49)  we  might 
distinguish  between  two  regions  of  the  location  of  the  interference 
signal's  direction; 

1)  |tan(0^/2)|  >  1/2  then  (49)  is  satisfied  for  any  interf erence- 
to-threshold  power  ratio;  however,  it, 

2)  |tan(6^/2)|  <  1/2  then  (49)  is  satisfied,  and  the  "accurately 
constrained"  canceler  performs  better  than  the  interference  canceler 
with  "ideally  filtered"  reference  loop,  only  if 


,  2-cos2(6./2)+/cos2(e./2)+cos4(8./2) 

<  _ i _ ± _ 1  _ 

2(5cos2(9i/2)-4] 


Now,  one  can  easily  show  that  for  -1/2  £  tan  (0^2)  5  l/2(l?.cos2(0i/2)l4/5) 

2 

the  left-hand  side  of  (50)  is  a  decreasing  function  of  cos  (6^/2). 
Therefore,  the  interference-to-threshold  power  must  be  bounded  above  by 
an  increasing  function  of  | 0^| .  In  particular,  for  6^  ■  0,  (50)  becomes 


8  1 1 1 2  1  1,207 


(51) 


In  Figure  6  we  sketch  the  different  regions  of  superiority  of  each 
canceler  as  a  function  of  8  . 

To  see  the  level  of  superiority  one  canceler  has  relative  to  the 
other,  we  write  from  (31)  and  (46) 


[  l+4g, I ■ “sin“(9^/2) 


F  /p  «  - ~ - i - - - - - 

d°  R°  cos2',ri/2)tl+4g;i:2(l+g;il2)sin2(0i/:)] 
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(Fdo/FRo}dB  =  201og[l+4g|l|2sin2(e./2)) 

-101og[l+4g|l|2a+g|l|2sin2(8i/2)]-20log[cos(ei/2)]  (53) 


Comparing  the  Desired  Signal-to-Noise  Ratio 
From  (32)  and  (A-31)  we  write 


1 S  |  2out  >  ^  1  S  1  “ 
Nout  2 


cos2 (9g/2) [ l+2g 1 1 j  2sin0^ (tan(6^/2)-tan(6g/2) ] 2 
l+4g[(|  s|  2(l-coses)+l  1 1  2(1-cos6  . ) )  (l+2g(  1 1|  2+|  S|  2)  ] 


(54) 


If  the  desired  signal  angle  6g  is  not  so  small,  and  the  relative  directions 
between  the  desired  signal  and  interference  as  well  as  their  powers  are  such 
that  the  second  terms  in  the  bracket  of  both  the  numerator  and  the  denominator 
of  (54)  are  much  less  than  one,  then 


1 S 1 2out  > 
M  out 


cos2(8  /2) 
aN 


(55) 


On  the  other  hand,  if  these  terms  in  the  brackets  are  much  greater  than  one. 


Chen 


i s  l.2  °u.c  1 2  iiii  . 

N  out  2 


!ll4sin2(ei/2)sin2(9i/2-6g/2) 

[  :i|2+|S|2][  !s|2sin2(5s/2)  +  |l!2sin2(6i/2)) 


( Sb'* 


Is'2 


j-  (l-cos(9i-9  ))  III  >>  ! S i ,  9g<9. 

N 


'  S '  2  !  I  j  ^  ,  ,  , 

— i—  - — lt-  (l-cos(6  ))jS|«  :  I ,  ,  9,  <9 

■  Ci  4  X  S  i  S 

"N  |S  1 


Except  for  the  factor  cos‘(“g/2)  in  (55),  these  are  the  same  results  as  the 
ones  obtained  with  the  interference  canceler  with  reference  signal  loop  when 
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it  operates  in  the  acquisition  mode  [3]. 

For  the  case  when  is  small,  and  in  particular  with  6s=0,  (32) 
becomes  (see  also  (A-30)) 


I  S ' 2out 
N  out 


[s|2  [l+4g!i;2sin2(ei/2)]2 

c2  l+8g|l!2(l+2g!l|2)sin2(e./2) 


If  2g | I | 2sin2 (e^/2) 


>>  1,  then 


(57) 


1  S  1  “out 
N  out 


2 


sin2(6./2) 

aN 


ill 


(1-cosS  . ) 


(58) 


Therefore,  when  the  direction  of  arrival  is  such  that  Icos6J<<1,  the  output 
SNR  is  approximately  equal  the  input  SNR  regardless  of  the  desired  signal 
power.  It  degrades  when  cose^  approaches  one.  This  result  is  the  same  as 
the  one  obtained  for  the  interference  canceler  with  reference  signal  loop 
when  it  operates  in  the  tracking  mode  [3]. 


CONCLUSION 

For  the  two-element  array  directionally  constrained  canceler  we 
establish  the  dependence  of  the  SIR  and  SNR  on  the  STR  and  the  ITR  power 
ratios  as  well  as  on  the  directions  of  arrival  of  the  desired  signal  and 
interference  relative  to  the  direction  of  constraint  (assumed  zero).  From 
this  we  obtained  the  improvement  factor  in  STR  and  SNR.  Two  modes  of 
operation  were  distinguished.  Firstly,  we  took  the  direction  of  arrival  of 
the  desired  signal  to  be  far  away  from  that  of  the  constraint  direction 
(6g  is  sufficiently  large)  and  found  the  improvement  factor  in  SIR  for  the 
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extreme  conditions  when  both  signal  powers  are  large  or  both  are  small. 
Secondly,  assuming  the  desired  signal's  direction  of  arrival  is  sufficiently 
close  to  the  constraint  direction  (8g  is  small),  we  establish  the  sensi¬ 
tivity  of  the  performance  of  this  canceler  to  the  inaccuracy  in  constraint 
(0  ^  0) .  We  found  that  for  sufficiently  small  9g  the  improvement  factor 

(in  dB)  is  quadratic  in  6g  with  a  maximum  at  a  point  which  might  differ 
from  9g  ■  0  (accurately  constrained).  Depending  on  the  direction  of  arrival 
of  the  interference  and  on  the  STR,  this  maximum  may  occur  when  the 

desired  signal's  direction,  9  ,  has  the  same  or  opposite  sign  (the  direc- 

sm 

of  constraint  is  taken  to  be  at  zero)  as  that  of  the  interference  direction. 

In  comparing  the  performance  of  the  "directionally  constrained  canceler 

with  that  of  the  "IMS  canceler,"  firstly,  in  relating  the  first  canceler 

when  9g  is  sufficiently  large  to  the  second  canceler  operating  in  the 

acquisition  mode,  it  was  shown  that  in  both  extreme  power  conditions 

(large  or  small  signals  powers),  the  first  canceler  is  superior.  Secondly, 

we  related  the  first  canceler  when  9  is  small  to  the  second  canceler 

s 

operating  in  the  tracking  mode.  In  particular,  special  attention  was  paid 
to  the  ideal  case  of  both  cancelers;  namely,  the  first  canceler  accurately 
constrained  and  the  second  canceler  having  an  ideally  filtered  reference 
loop.  It  was  found  that  the  directionally  constrained  canceler  performs 
better  with  any  ITR  only  if  the  interference  direction  of  arrival  is  such 
that  ]tan9^/2  >  1/2  (of f -main-beam  interference).  If,  however,  the 
interference  gets  closer  (in  main-beam  interference)  then  the  direction- 
ally  constrained  canceler  is  preferred  only  if  the  ITR  is  smaller  than  a 
certain  value.  This  value  becomes  smaller  and  smaller  when  the  interfer¬ 
ence  direction  of  arrival  gets  closer  and  closer  to  the  constraint  direction. 
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This  behavior  suggests  complementary  properties  of  the  two  cancelers, 
so  that  still  better  performance  can  be  obtained  if  we  use  both  directional 
(directionally  constrained)  and  spatial  (reference  loop)  information  of  the 
desired  signal.  The  SNR  performances  of  both  cancelers  are  found  to  be 
quite  similar. 
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c-i 

Using  (17)  and  (1)  we  have 

Y*(t)Y1(t)  -  |tS1(t)-S2(t))+(I1(t)-I2(t))+(N1(t)-N2(t))J l2 


-  ]S1(t)-S2(t)l2  +  | I1(t)-I2(t) 1 2  +  |N1(t)-N2(t)|2 

-  Is^i2  +  +  2o2  a-u 

where  from  (3)  and  (4)  we  assumed  Chat  js(t)|  •  1,  |i(t)|  “  1;  and  that  the 

signals  and  noise  processes  are  mutually  uncorrelated. 

Similarly,  from  (16)  and  (17) 


Y^OYqU) 


[(S1(t)-S2(t))  +  (Il(t)-I2(t)  +  N1(t)-N,(t))]* 


•  {(S1(C)+S2(C>)  +  (Il(t)+I2(t))  +  (N1(t)+N2(t))] 


(s1-s2)*(s1+s2)  +  (i1-i2)*(i1+i2) 


The  noise  terms  cancel  out  since  we  assumed 


2 

'N  ' 


(C-2) 
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Substituting  for  Y^(t)  from  (16)  together  with  (1),  and  for  W  from  (19) 
we  obtain 


f:- ;>  i 


VC)WV 


/ 


-g[S1-S2)*(S1+S2)+(I1-I2)*(I1-M.2)] 

i+g(|s1-s2i2+;i1-i2;2) 

[(S1(t)-S2(t))  +  (N1(t)-N2(e))+  (I1(t)-I2(t))j 

-g[  |s1-s2!2(s1(c)+s2(c))+(i1-i2)*d1+i2)(s1(t)-s2(t))] 
i+g(is1-s2l2+li1-i2|2 


-g[  !i1-i2|2d1(t)+i2(t))+(s1-s2)*(s1+s2)d1(t)-i2(t)] 
i+g(js1-s2i2+|i1-i2'2 

-g(  (s1-s2)*(s1+s2)+d1-i2)*d1+i2)](N1(t)-N2(t)) 
i+g(is1-s2f2+|i1-i2;2 


Together  with  (17)  and  (1)  we  have,  after  simple  algebraic  manipulation 


Yo(t)+Y1(t)Wy 


(i+g!i1-i2|2)(s1(t)+s2(t))-gd1-i2)*d1+i2)(s1(t)-s2(t)) 
l-t-g  c ;  s1-s2 1 2+ 1  i1-i2  ! 2) 


+  — 


(l+g|S1-S2|2)(I1(t)+I2(t))-g(S1-S2)  (S1+S2)(I1(t)-I2(t)) 


i+g(|s1-s2!2+!i1-i2;“) 


+  N1(t)+N2(t)- 


g[(S1-S2  )(S1+S2)+(I1-I2)  (I1+I2)](N1(t)-N2(t)) 

i+g(!s1-s2!2+|i1-i2!2) 


S1(t)[l+2g(,lI2!2-I*I2)l+S2(t)[l+2g(|l1'|2-I2I1)] 

i+g(|s1-s2j2+U1-i2!2) 

I1(t)[l+2g(|s2j2-s*s2)]+s2(t)[l+2g(!si:2-s2si)] 

i+g(|s1-s2!2+!i1-i2|2) 

N1(t)[l+2g(|S2|2-S*S2+!l2i2-I*I2)]+N2(t)[l-f2g(,S1;2-S*S1-t-  ^ 

i+g(!s1-s2i2+|i1-i2'2)  (C_4) 
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Using  (3)  and  (23) 


j6s,  ^s. 


S1(t)+S2(t)  =  s(t) | S | (e  1  +  e  2) 

2(6s  +9s  } 

-  2s(c) | S !  e  1  2  cos(6s/2) 

j6  j (6 .  -0  ) 

I1*I2S1(t)  «  s(t) | S ( e  1 1 1 1 2 e  2  1 

1  6s 

f(ec  +e  )  i(-y  -  e  ) 

7  ^  s?  ^  1 

-  s(t)|S| |l|e  2  e 


where  we  also  used  (24).  Also 


I2I1S2(t)  -  s(t) | S | | I |  e 


4(9  +9  )  j(6  -  -f) 

2  2  si  s2  i  2 


Hence , 


4(9g  +6^  ) 

I*I2S1(c)+I*I1S2(t)  -  2s(t)]s| |l|2e  1  32  cos(0.-es/2) 


Similarly,  2(Si  +9i  } 

l1(t)+I2(t)  -  2i(t)|l| |s|2e  1  2  cos(91/2) 

2(9i  +9i  > 

S*S2I1(C)+S*S1I2(C)  -  2i(t)|l| | sj  2e  1  2  cos(es-e1/2) 


Also 


1 2  * 

Sl  -S1S2  ' 

~  -j  9 

| S | 2(l-e  S) 

(C-9 ) 

1 2  * 

,  ,2,  -j9i. 

Xl  * 

| I l  (1-e  x) 

(C-10) 

IS1-S2|2  " 

2  j  S 1 2 (1-cosfi  ) 
s 

fC-11) 

! 1i_I2l 2  " 

2 ! X 1 2(l-cosei) 
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Equation  (38)  can  be  written  in  the  form 


F.  =  F,  •  L. 
d  do  d 


(C-13) 


where  represents  the  change  in  the  improvement  factor  due  to  a  non-zero 
6g  (i.e.,  non-accurate  constraint).  This  change  in  improvement  factor  can  be 
written  (in  dB)  as 

(I'd) dB  *  201ogcos(@g/2)+201og(l-a2)-201og(l+a2)  (C-14) 

where 

tan(8  /2) 

al  *  tan(e./2)  “  (c_15) 


a,  -  4g|s|2sin2(6  /2)(l-tan(e  . /2) /tanO  12))  (C-16) 

^  s  i  s 

For  i 6 s |  <  |0  |  |a^|  <  1.  Also,  one  can  easily  show  that  for  ja2|  <  1  it  is 

sufficient  to  require  that  0g  be  small  enough  to  satisfy 


sin2  (8  /2)  +  sine  <  l/4g  |  S  |  2 

s  2i/2"  s 


( C— 1 7  ) 


where  we  took  |0^|  <  tt/2. 

Under  these  conditions  we  can  use  the  series  representation  of  log(l-a^), 
log (l+a2),  and  log  cos(6s/2)  and  write 


(VdB  =  -20[’ 


(0/2)' 


(9  nr 
+  —4^ —  + 


-20(a2-  -f  + 


tane  H  fane  / 2 

-20  [ - 2_  +  - 2 —  4 

tan8./2  2  . 

1  t  anc  .  2 
1 

(C-18) 


To  a  first  approximation  the  change  in  improvement  as  a  result  of  a  non- 
accurate  constraint  is  given  by 


024 


ssaaa 


i 


r 


,6  /2  2  tan6  /2  o  i 

(Ld)dB  =  "20[  “2~  +  +  ^elslsin  (8g/2)(l-tan(ei/2)/tan(8s/2))] 

Since  9  is  small,  sin(9  /2)  =  tan(9  /2)  *  9  /2,  (A-18)  can  be  written  as 
s  s  s  s 

[Ld]dB/-10  =  (es/2)2[l+8g|si2]+2(es/2)[l/tan(ei/2)-Ag|s|2tan(e./2)] 
The  left-hand  term  of  (A-19)  has  a  minimum  at 


9  /2 

sm 


[l/tan(ei/2)-4g| S|  can(ei/2) ] 

l+8g|s!2 


The  value  of  this  minimum  is 


tl/tan(9i/2)-4g|  S|  2 tan (9^2)  ]2 


l+8g|s|' 
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Using  trigonometric  relations  we  write  (49)  as 


[l+2g | l| 2(l-cos0  ) ]2  ? 

- —  >  l+2g  1 1 1  2  ( 1+g  l!  )(l-cos9.) 


l+4g | I | 2(l-cos9  )+4g2 | I | 4(1-cos9  ) 2  . 

>  - __i  [i+2gjl!2(l+g|lj2)(l-cos9.)] 


Arranging  terms  we  get 

2 

«  ,  ~  1-cos  9.  ~  1+cose . 

g  ji:  (1-cosS  )  [4 - i — 5-]-^2g  1 1 !  (l-cos9  . )  [  2 - 5 — — ] 

1  (l-cos8.r  1 

+  -|(l-cos  )  >  0 


2 

ax  bx  +  c  >  0 


where  . 

x  *  g ! I i  (1-cosS^)  i  0 


(C-18) 

(c-19; 

(C-20) 

( C-21 ) 


(C-22) 


(C-23) 

(C-24 ) 


r:-2S 


l+cos6 .  , 

a  =■  4  -  — - —  =  4  -  cotn"(9./2) 

l-cos6 .  1 


(C-25 ) 


l+cos0 . 


b  =  2[ 2- 


2(2-cos  (6  /2)) 


( C— 2  6 ) 


c  =  j  (l-cosei)  =  sin2 (9^/2) 


( C— 2  7 ) 


Thus, 


2  6  i. 


(•j)2  -  ac  «  (2-cos2(6  ./2))2-sin2(6  ,/2)[4-cotn"  ~ 


cos^*(8i/2)  +  cos2(ei/2) 


(C-28) 


We  are  looking  for  a  positive  solution  x  of  (A-23).  Notice  that  c  £  0; 
therefore , 

1)  if  a  >  0  or  |tan(9^/2);  >  1/2,  then  every  x  >  0.  (c-23) 

is  satisfied,  however, 


2)  if  a  <  0  or  tan  (%J2)  <  1/2,  then  it  is  easy  to  show  that 


onlv  when 


x  »  g 1 1 j  2(l-cos6i)  s 


-(2-cos2 ( 6 i/2>  —  / :os^(9^/2)+cos2(Si/2) 
4-cotn2 (6i/2) 


(C-23)  is  satisfied.  Or,  equivalently,  if 


8 !  I 


2  2-cos2(8i/2)+  /cos2(9i/2)+cos4(6i/2) 

<  - 

2(5cos^(ei/2)  -  4) 


(C-29) 
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From  (32) 

| l+2g i S | 2 (1-e  S)+| X | 2(l-e  i))|2 

-  l+4g(!s|2(l-cosSg)+|l|2(l-cos9i))+8g2[ |s|2!l|2(l+cos(eg-9.)-cose.-coses) 
+  [ s( 4(l-coseg)+[ 1 [ 4(l-cosei) ]  (C-30) 
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l+4g(|s|2(l-cos6s)+| lj2(l-cosei) 

|s|2|l|2(l+cos(5  -8 ,)-cos6 ,-cosS  )+| S | 4(l-cos6  )+|l|4(l- 
[l+2g  - i— ^ 1 5 - 

| S |  (l-cos8g)+| I |  (l-cos6i) 

l+4g ( ! S I  2 ( 1-C0S9  s)+ 1 1 1  2 (1-cose ±) ) [ l+2g ( I S I 2+ 1 1 ( 2 ]  (c-31 ) 
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Figure  3  Directionally  Constrained  Interference  Canceler 
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gure  5  Two-element  Interference  ConcelerWith  Directional  Constraint 
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EFFECT  OF  POINTING  ERROR  ON  A  DIRECTIONALLY-CONSTRAINED  ADAPTIVE  ARRAY 

We  are  here  concerned  with  a  satellite-borne  adaptive  array  for  the 
uplink  of  a  communication  satellite  serving  widely  scattered  multiple  users. 
Ground  positions  will  be  known,  at  least  approximately,  making  it  reasonable 
to  consider  a  directionally  constrained  adaptive  algorithm  such  as  that 
described  by  Applebaum  and  Chapman  [1].  Of  concern,  however,  is  the  effect 
of  pointing  error,  which  may  arise  because  of  imperfect  knowledge  of  ground 
position,  or  because  of  error  in  determining  satellite  orientation,  or  because 
of  diffraction  effects.  An  investigation  of  the  sensitivity  of  the  signal-to- 
interference-plus-noise  ratio  (SINR)  to  error  of  this  kind  in  the  case  of  the 
Applebaum  array  with  zero-order  directional  constraint,  was  carried  out  and 
is  here  reported. 

The  analysis  was  based  on  a  linear  array  with  randomly  deployed  elements. 
For  numerical  results  a  10-wavelength  array  was  used,  both  filled  and 
sparse.  The  array  structure  apt  to  be  used  in  the  final  application  is  likely 
to  be  different  both  in  form  and  size,  but  the  phenomena  under  study  is 
expected  to  follow  a  similar  pattern. 

The  array  processor  is  shown  in  one  of  its  forms  in  Figure  1.  If  the 
array  were  linear,  ^s  shown  in  Figure  2,  and  if  a  desired  signal  were  arriving 
from  direction  6^, the  induced  voltage  in  the  ith  element  would  be  the  real 
part  of 

j lw  (t+T  )  +  $] 

vi(t)  =  a(t+T1)e  C  (1) 

where  =  x^cos8j/c,  c  is  the  propagation  speed  along  the  ray  path,  a(t) 

(1]  S.P.  Applebaum  and  D.J.  Chapman,  "Adaptive  Array  with  Fain  Beam 

Constraints,"  IEEE  Trans.  Anntenas  and  Tropag.,  Vol.  AP-26,  September 

1976. 
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is  the  complex  baseband  signal,  ujc  is  the  radian  carrier  frequency,  $  is  the 
carrier  phase  at  the  reference  point,  and  is  the  position  of  the  ith 
element  relative  to  a  reference  point  (origin). 

The  signal  at  the  reference  point  itself  is  taken  to  be 


v  (t) 
o 


j  (uj  t  •+  <}>) 

a(t)e 


(2) 


In  this  work  the  baseband  signal  a (t)  is  assumed  to  be  of  narrow  bandwidth, 
changing  only  slightly  in  the  interval  T  for  all  i  =  1,2,...N.  Thus,  we 
assume  aU+t^)  =  a(t)  and 


jkx  cos0 

vt(t)  =  a(t)  e 


)  (u^t-hfi) 
e 


j  (u)  t-H}>) 

“  vt(t) e 


,  i 


1,2, ...N 


(3) 


where  k  *  2ir/X  =  id  /c. 

c 

In  the  following  we  work  explicitly  with  the  complex  coefficient  v  ^.'c)  , 
suppressing  the  carrier  factor  (w£t+<!>)  ,  Furthermore,  these  N  complex 

quantifiers  are  represented  by  the  vector  v_  -  (t)  »v  >  *  *  -VN  (t))  .  Noise 

and  interference  are  similarly  treated. 

The  beamformer  part  of  the  processes  Is  comprised  of  a  set  of  phase 
shifters  and  an  adder  which  acts  to  form  e  -  j>  v  where  S  is  a  pointing 
vector  (T*  signifies  transpose  conjugate)  given  by 


jkx^  cose 
e 

JkXjCosS 


jkxN 


cose 


(4) 


e  is  the  beam  pointing  direction.  If  9  -  6d,  the  direction  of  the  desired 

signal  arrival  e  *  Na(t). 

c& 

Simultaneously,  the  phase  shifted  signals  are  applied  through  an 
Nx(N’-l)  matrix  transformer  A,  chosen  such  that  A£  “0,  to  a  multiple  side- 
lobe  canceler.  The  matrix  transformer  forms  a  set  of  (N-l)  outputs  which 
contain  no  signal  from  the  pointing  direction.  If  the  signal  magnitude  is 
the  same  at  each  element,  this  device  only  needs  to  form 


visl  "  vi+lsi+l 


1=1,  2. . . (N-l) 


(5) 


The  multiple  sidelobe  canceler  (see  Applebaum  [1],  Figure  1)  acts  to  form  can¬ 
celling  signals  to  those  inputs  arriving  from  directions  other  than  the  point 
ing  direction  9.  Applebaum  points  out  that  the  scheme  approaches  the  maximum 
SlNR  combiner  defined  by 


v 

o 


(6) 


where 


W 


-1  * 

M  AS 

T  -1  "* 
S1!*  S 


(7) 


and  M  is  the  covariance  matrix  of  the  complex  signals  arriving  at  each  element 
that  is, 


*  T 

M  =  <  v  v  > 

The  SINR  at  the  array  processor  output  is 


(8) 


SINR  -<|a2 |>£T  M-1  S*  (9) 

2 

where  |a  |  is  the  mean  square  power  at  each  sensor,  provided  the  desired 

signal  arrives  from  the  pointing  direction. 

Now,  if  the  pointing  directiot  is  different  from  the  arrival  direction 

of  the  desired  signal  the  SINR  wi'I  be  less  than  that  given  by  (9).  An 

expression  for  the  SINR  in  this  case  is  now  obtained.  With  9“0,,  the  desired 

th  ° 

signal  arrival  direction,  and  6  •  ^  undesired  signal  direction 

we  define  direction  vectors  and  S^,  Also  defining  a  noise  vector  n  of 
independent  components,  we  have 


M  -  <|([aSd  + 


I.ek 
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where 
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<JB^|>“  mean  interference  power  of  the  kC^  interferer  at  each  sensor 
2 

0  *  noise  power  at  each  sensor 

and  the  8^  for  the  K  interferers  are  independent  complex  random  processes. 

0  and  the  corresponding  pointing  vector  £  (without  subscript)  will 
refer  to  an  arbitrary  pointing  direction  and 

-1  * 

W  =  C  M  S  (restatement  of  (7)) 

will  represent  the  equivalent  weights  generated  by  the  processor  when  point- 

T  — 1  if 

ing  in  the  direction  9.  C=(SM  S)  is  a  constant  for  a  given  pointing 

angle.  The  output  power  of  a  desired  signal  arriving  from  direction  0.  is 

d 

pd  -  <(WT  a  Sj)*  (a  sj  W)>  -  <a2  >C2  sT  M-1  Md  M-1  S* 

Similarly  the  output  interference  power  of  the  kC^  interferer  is 

?Ik  -  <Sk>  °2  A"  M'1  M*  M-1  S*  (12) 

and  the  output  noise  power  is 

2  2  T  -1  * 

PN  -  a <T  S1  M  A  IM  A  S  (13) 

the  output  SINR  Is 


2  T  -1  —1  * 

<a  >  S  M  M .  M  S 

—  d  — 


SINR  - 


J1<6k  >  -  M"1  Mik  M_1  £*  +  °2  ST(M-1)2  S* 


Equation  (1A)  has  been  used  to  calculate  the  SINR  for  a  number  of  specific 
cases  as  described  below. 


(a)  Filled  Array. 

An  equally- spaced- element  linear  array  was  assumed;  the  interspace  between 
elements  is  y  and  the  beamwidth  of  this  array  is  approximately  y  m  ^  rad  -  5.7* 
The  desired  signal  is  assumed  to  arrive  from  9^  -  90°  and  a  single  inter¬ 
ference  signal  is  assumed  to  arrive  from  0^  ■  85°.  The  signal  power  is  chosen 
to  be  10  dB.  The  noise  power  is  0  dB  ,  and  the  interference  power  is  chosen 
as  10  db.  5  db.  20  db.  By  varying  the  pointing  angle  0  (the  angle  to  which  the 


! 


array  is  focused),  the  results  shown  in  Figures  3,  A,  5  are  obtained.  At  9 
■  90°,  meaning  that  the  array  is  focused  exactly  toward  the  arriving  desired  sig¬ 
nal,  the  values  of  SINR  for  the  three  levels  of  interference  are  about  equal.  This 
jneans  the  array  nulling  is  proportional  to  the  interference  power.  In  fact,  the 
curves  of  Figures  3,  A,  5  are  about  the  same  except  near  6  =  85®,  the 
direction  of  interference  arrival.  When  the  array  is  pointed  at  this  angle  the  in¬ 
terference  power  is  constrained  to  be  constant  while  the  signal  power  is  minimized. 
(SEE  NOTE  BELOW) 

(b)  Sparse  Array. 

A  seven-element  sparse  array  of  length  10X  was  simulated  as  follows.  A  uniform¬ 
ly  distributed  random  number  generator  was  used  to  select  seven  numbers.  These  num¬ 
bers  were  then  linearly  scaled  to  fill  the  range  (-5X,  5X)  so  that  one  element 
would  appear  at  x  *-5X  and  another  element  would  appear  at  x  =5X.  One  such  array 
was  simulated  with  positions  which  turned  out  to  be  -5. OX,  -3.A7X,  -1.55X,  -1.17X, 
0.2AX,  A.33X,  5. OX.  The  beanwidth  of  this  thinned  array  is  about  the  same  as  that 
of  the  21-element  equal-spaced  array.  Numerical  results  were  obtained  as  described 
below.  Signal  is  always  assumed  to  arrive  from  9d“90°.  (SEE  NOTE  BELOW) . 

I.  One  interferer  at  Gj-  85°  gave  the  result  shown  in  Figure  6.  Compared 

with  Figure  3.  the  SINR  is  5dB  less  at  9-90°,  but  it  is  not  so  sensitive 
in  the  neighborhood  of  9-90°.  The  loss  in  gain  is  a  direct  result  of  the 
use  of  fewer  elements  here  than  in  the  case  of  the  filled  array. 

Ii.  Figure  7  is  the  result  of  two  interferers  at  0jj*85O  and  6^^95°.  The 

SINR  at  90°  is  about  the  same  as  that  in  Figure  6,  but  it  is  less  sensi¬ 
tive  around  90°  than  Figure  6. 

III.  To  see  the  effect  of  the  nulling  in  the  main  beam,  we  assumed  one  interfer¬ 
er  coming  in  from  9^89°.  The  result  is  shown  in  Figure  8.  At  90°  the  SINR 
is  lOdB,  a  loss  of  about  8dB  in  SINR  with  no  pointing  error,  but  the  array 
is  much  less  sensitive  than  in  Figures  3  -  7.  We  then  assumed  that  there 

are  four  interferers  (at  85°  89°  91°  95°).  Two 'of  them  (89°  91°)  are  in  the 
main  beam  when  pointing  angle  is  90°.  The  result  (Figure  9)  shows 

NOTE:  Figures  3  through  18,  following,  are  plots  of  S/N+I  as  a  function  of 
pointing  angle.  Figures  3  through  5  are  derived  from  a  simulated,  21-element- 
equal-spaced,  linear,  filled  array  of  length  10X.  Figures  6  through  18  are  de¬ 
rived  from  a  simulated,  7-element-non-unif ormly-spaced ,  sparse  array  of  length  10X. 
Opposite  the  figures,  only  the  respective  measurement  conditions  arc  indicated. 

In  all  figures,  the  desired  signal  is  assumed  to  arrive  from  90°,  the  signal's 
power  is  lOdB,  and  noise  power  is  0UB. 


FIGURE  3.  1  INTERFERER,  POWER  lOdB, 
AT  85° 


FIGURE  4.  I  INTERFERER,  POWER  5dB, 
AT  85° 


FIGURE  5.  1  INTERFERER,  POWER  20dB, 
AT  85° 


FIGURE  6.  1  INTERFERER,  POWER  lOdB, 
AT  85° 


FIGURE  7.  2  INTERFERERS,  POWER  lOdB  EACH 
AT  85°,  95° 


FIGURE  8.  1  INTERFERER,  POWER  lOdB, 
AT  89° 


FIGURE  9.  4  INTERFERERS,  POWER  lOdB  ' 
AT  85°,  89°,  91°,  95° 


FIGURE  10.  4  INTERFERERS,  POWER  lOdB  EACH 
AT  37°,  127°,  190°,  253° 


yh+i  «b) 


FIGURE  11  4  INTERFERERS ,  POWER  20dB 

EACH,  AT  37°,  127°,  190°,  253 


FIGURE  11a  5  INTERFERERS,  POWER  lOdB 

EACH,  AT  37°,  65°,  127°,  190°,  253° 


FIGURE  11b  5  INTERFERERS,  POWER  20dB 

EACH,  AT  37°,  65°,  127°,  190°,  253° 


FIGURE  12  6  INTERFERERS,  POWER  lOdB 

EACH,  AT  37°,  65°,  127°,  190  ,  253°,  335° 


FIGURE  13  6  INTERFERERS,  POWER  20dB 

EACH,  AT  37°,  65°,  127°,  190  ,  253°,  335° 


FIGURE  14  7  INTERFERERS,  POWER  lOdB 

EACH,  AT  0,37,65,127,190,253,335  (DEG) 


FIGURE  15  7  INTERFERERS,  POWER  lOdB 

EACH, AT  37,65,127 ,190, 253, 268. 5, 335 (DEG) 


FIGURE  16  7  INTERFERERS,  POWER  lOdB 

EACH, AT  37, 65, 127, 190, 253, 268. 5, 335 (DEG) 


FIGURE  17  5  INTERFERERS,  POWER  lOdB 

EACH,  AT  65  ,  75°,  85°,  95°,  105° 


FIGURE  18  6  INTERFERERS,  POWER  lOdB 

EACH,  AT  65°,  76°,  85°,  95°,  105  ,  115° 


D-9 


POINTING  ANGLE  (DEGREES) 


the  array  is  not  able  to  handle  this  case.  The  sensitivity  to  pointing 
error  is  small  but  the  SINR  has  fallen  below  OdB  even  with  no  pointing 
error. 

IV.  We  arbitrarily  picked  four  angles,  37“,  127°,  190°,  253°  as  interfe¬ 
rence  arrival  angles  to  simulate  widely  spaced  interferers,  not  in  the 
main  bean..  The  results  are  shown  in  Figures  10  and  11.  In  Figure 

10  the  signal  power  is  lOdB.  Interference  powers  are  also  lOdB.  In  Fig¬ 

ure  11  signal  power  is  lOdB  and  interference  powers  are  20dB.  The 
results  as  indicated  by  Figures  10  and  11  are  about  the  same.  The 
SINR  at  90°  is  16.7dB  and  the  difference  between  this  and  the  result  in 
Figure  6  is  less  than  2dB. 

V.  Five  interferers  at  37°,  65°,  127°,  190°,  253°  were  simulated.  The  re¬ 
sults  art  in  Figures  11a  and  lib.  Again,  lOdB  and  20dB  interference 
powers  were  used.  The  results  are  about  the  same,  and  the  SINR  at  90° 
is  very  c:ose  to  that  in  Figures  10  and  11. 

VI.  Six  interferers  were  then  assumed  at  37°,  65°,  127°,  190°,  253°,  336° 
with  results  shown  in  Figures  12  and  13  for  interference  powers  of 
10  and  2 OdB  respectively.  The  results  are  much  different  from  those  of 
four  and  five  interferers.  The  three  major  differences  are:  (1)  the 
SINR  at  90°  is  about  6dB  less,  (2)  the  pointing  error  sensitivity  is 
less,  (3)  the  curves  for  different  interference  powers  are  quite  differ¬ 
ent.  For  most  angles  the  SINR  for  20dB  interference  power  is  larger  than 
the  SINR  for  lOdB  interference  power.  The  following  may  explain  these 
phenomena.  Since  the  array  consists  of  seven  elements,  it  can  null  six 
Inputs.  For  the  six-interferer  case,  if  we  point  the  array  at  an  angle 
from  which  no  signal  or  interference  comes,  the  array  will  try  to  null 
all  the  signal  and  interferers.  However,  the  array  can  handle  at  most 
six  inputs  and  there  are  seven  signal  and  Interferers  together.  The 
larger  the  interference  powers  are,  the  less  attention  the  array  will 
pay  to  the  signal.  It  is  similar  to  the  power  inversion  array. 

VII.  For  the  six-interferer  situation,  if  we  point  the  array  exactly  at  90° 
it  only  needs  to  null  these  six  interferers,  and  should  be  able  to 
handle  that.  If  there  are  seven  interferers,  when  pointing  at  90°, 
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the  array  is  not  going  to  work  well,  and  we  expect  the  SINR  at  90®  to 

• 

be  poor.  To  see  this  we  added  another  interferer  coming  from  0j  *=  0 
to  the  six  interferers  in  VI.  However,  the  simulation  result  shows 
that  the  SINR  at  90®  for  these  seven  interferers  is  onlv  1  db  less  than 
that  of  six  interferers  (compare  Figures  14  and  12).  .  In  order 
to  explain  this  phenomenon,  we  calculated  the  power  pattern  for  some 
of  the  cases  treated  above, believing  that  the  array  with  6  interferers 

■ay  have  a  natural  null  at  0*  when  it  is  focused  at  0  -  90®.  The 

T  2  T 

power  pattern  of  the  array  is  P(B)  -  |w  sj  ,  where  W  is  the  weight 
vector  when  signal  arrives  from  90®  and  S  is  the  steering  vector 
which  we  vary  from  0®  to  360®.  By  the  constraint  we  have  P(90®)  »  1. 
The  power  patterns  are  given  in  Table  1. 


Direction 

5  Interferences 
(37®,  65®, 

127®,  190”,  253®) 

6  Interferences 
(37®,  65®,  127® 
190®,  253®.  335°) 

7  Interferences 
(0.  37®,  65®,  127 
190”,  253®,  335® 

0® 

*  -8.29  dB 

*  -9.97 

dB 

-23.33  dB 

37® 

-38.7  dB 

-28.40 

dB 

-32.34 

65® 

-46.12 

-24 . 64 

-24.63 

127® 

-41.55 

-19.69 

-18.30 

190’ 

-39.88 

-27.78 

-26.14 

253 

-43.75 

-27.43 

-21.20 

333 

*  -6.21 

-18.45 

-15.90 

TABLE  1  Power  Pattern  when  Weights  are  Fixed  for  Signal  Arriving  from 

0.  ■  90®,  and  5,  6,  and  7  Interferers  Arrive  as  Specified. 

d 

The  patterns  above  are  based  on  signal  and  interference  powers  of  10  dB 
each,  and  noise  power  is  0  dB.  The  numbers  with  an  asterisk  (*)  in  front  mean 
that  'ihere  is  no  interference  coming  from  those  directions. 

From  Table  5.1  we  can  see  that  even  for  6  interferers,  the  array  has  a  fair¬ 
ly  good  null  at  0®,  so  that  when  another  interferer  arrived  from  0®,  the  array 
naturally  suppressed  it.  We  therefore  chose  two  other  angles,  280®  and  268.5®, 
where  the  values  of  the  power  pattern  with  6  interferers  at  these  angles  are 
-3.91  dB  and  2.94  dB  respectively.  With  the  seventh  interferer  coming  in  from 
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these  angles,  the  SINR  at  90°  is  8.82dB,  or  3.75dB.  Thus,  the  array  does  not 
handle  these  cases  as  well  as  the  one  in  which  the  seventh  interferer  arrives 
from  0  .  But  the  results  are  still  not  decisive.  We  may  argue  that  in  the  case 
of  seven  .or  even  more  interferers  the  array  processor  still  seeks  to  minimize 
the  array  output  to  all  inputs  except  for  the  one  in  the  constraint  direction.  If 
it  can  only  null  in  a  smaller  number  of  directions  it  places  them  in  such  a  way 
that  the  larger  number  of  unwanted  inputs  is  minimized*,  thus  the  result  observed 
for  7  interferers  should  not  be  surprising.  What  is  perhaps  less  obvious  is  why 
the  SINR  is  not  better  with  6  interferers.  For  one  to  five  interferers,  the  SINR 
when  pointing  toward  the  desired  signal  is  always  within  about  2dB  of  the  best  ob¬ 
tainable.  For  the  6-interferer  case,  it  falls  8  to  lOdB  lower.  However,  this  may 
be  a  fortuitous  result  of  the  particular  choice  of  array  element  layout;  the  next 
article,  where  we  try  different  layouts,  shows  that  interference-cancelling  for 
some  layouts  may  be  poorer  than  for  others. 

VIII.  In  III  we  simulated  a  case  of  main  beau.  nulling,  and  it  showed  that  the 
array  did  not  work.  In  IV-VII,  we  the  efore  spread  out  the  interferers 
and  saw  acceptable  results.  It  is  interesting  to  see  what  happens  when 
multiple  interferers  are  compact  but  still  not  in  the  main  beam.  There¬ 
fore  we  chose  interference  arrival  angles  of  65°,  75°,  85°,  95°,  105° 
for  5  interferers  and  65°,  75°,  85°,  95°,  105°,  115°  for  6  interferers. 

The  results  plotted  in  Figures  17  and  18  are  about  the  same.  The 
SINR  is  at  least  as  good  at  90°  as  for  the  scattered  interferer  case  (it 
is  much  better  for  6  interferers)  and  the  sensitivity  to  pointing  error 
is  better. 

The  results  of  this  analysis  and  computation  indicate  that  pointing  sensiti¬ 
vity  is  such  that  pointing  error  of  the  order  of  0.1°  may  significantly  reduce 
the  output  SINR.  Accurate  control  of  pointing  is  therefore  viewed  as  being  essen¬ 
tial.  It  also  confirms  the  importance  of  making  the  array  large  enough  to  put  po¬ 
tential  interferers  outside  the  beamvidth  of  the  array  beamformer  in.  Figure  1. 
Very  likely,  the  sensitivity  to  pointing  error  will  Increase  with  increasing  array 
size,  making  it  even  more  necessary  to  incorporate  some  mechanism  for  accurate 
pointing. 
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In  an  earlier  section  of  this  report  ("Effect  of  Pointing  Error 
on  a  Directionally-Constrained  Adaptive  Array,  Appendix  D),  computational 
work  reported  on  SINR  sensitivity  to  pointing  error  was  based  on  a 
single  random  deployment  of  array  elements.  As  part  of  our  study 
of  the  effect  of  array  structure  we  have  continued  the  investigation 
reported  above,  determining  the  SINR  for  various  random  element  deployments. 

The  array  was  linear  and  comprised  of  seven  elements;  the  powers 
of  signal,  a  single  interferer,  and  noise  were  taken  to  be  lOdB,  lOdB, 
and  OdB,  respectively.  Signal  was  assumed  to  arrive  broadside  (Bj  =  90°) 
and  interference  was  assumed  to  arrive  from  various  angles  0j.  The  pointing 
vector  was  set  to  look  in  the  direction  of  the  signal. 

Fifty  different  layouts  were  selected  using  a  random  number  generator.'  Two 
different  cases  were  treated.  In  the  first,  element  locations  were  based  purely 
on  a  uniform  distribution  over  (-_5X,  5X).  In  the  second  case  the  same  random  dis¬ 
persion  of  elements  was  scaled  to  fully  cover  the  interval  (-5X.  5  X)  •  The  results 
are  presented  in  Figures  1  and  2  for  unsealed  and  scaled  arrays,  respective¬ 
ly.  The  A  situations  shown  in  Figure  1  are  for  Interference  arriving  from  di¬ 
rections  10°,  5°,  2°,  and  0.1°  from  the  signal  direction.  The  first  two  simulate 
interferers  at  or  beyond  the  beamwldth  of  the  largest  arrays  in  the  50  samples, 
the  other  two  simulate  interference  generally  inside  the  beamwidth.  We  draw  the 
following  conclusions  from  these  results.  For  interference  adequately  spaced  away 
from  the  desired  signal  the  values  of  SINR  are  largely  concentrated  near  18.45dB, 
the  maximum  SINR  obtainable.  Though  some  cases  are  found  in  which  the  SINR  is 
poorer  by  about  3dB,  these  cases  are  rare.  The  loss  in  SINR  must  be  associated 
with  imperfect  nulling  arising  perhaps  from  poor  element  placement.  When  inter¬ 
ference  is  within  the  beamwidth  of  the  array,  the  values  of  SINR  are  more  uni¬ 
formly  spread  and  almost  always  lower  than  it  is  for  interference  outside  the 
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beamwidth.  The  constraint  on  the  main  beam  gain  makes  it  virtually  impossible  to 
get  a  good  null  on  the  interference.  One  should  note  that  interference  about  half¬ 
way  into  the  beamwidth  (2°  from  the  desired  signal)  produces  a  median  S1NR  about 
5dB  poorer  than  the  best  obtainable.  For  interference  very  close  to  the  desired 
signal  (the  0.1°  case)  the  SINR  hovers  around  0.2dB,or  more  than  18dB  poorer  than 
the  best  obtainable. 

Because  some  sample  arrays  will  not  span  10X  when  element  positions  are 
randomly  chosen, we  scaled  the  arrays,  stretching  them  to  fill  the  10A  space. 

Two  cases  of  interference  arrival  were  simulated,  10°  and  2°  from  the  desired 
signal.  For  the  10°  case  the  values  of  SINR  generally  tended  to  be  nearer  the 
best  obtainable  value  but  there  still  were  cases  about  3dB  poorer.  The  2°  case 
was  not  significantly  different  than  the  unsealed  case. 
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Adaptive  interference  cancelling  arrays  utilizing  either  known 
direction  of  signal  arrival  or  known  signal  waveform  structure  have  been 
analyzed  extensively  in  the  past  [1,2] .  The  former,  called  here  the 
directionally  constrained  array,  is  the  more  natural  one  to  use  in  point- 
to-point  communication  though  the  sensitivity  to  pointing  error  is 
substantial  [3,4].  Variants  of  the  latter  canceller,  in  which  the 
waveform  reference  is  generated  in  the  receiver  using  some  known  signal 
characteristic  such  as  the  code  used  in  spread  spectrum  modulation  has 
also  been  advanced  for  communications  applications.  We  refer  to  these  as 
IMS  interference  canceler  arrays. 

Recently,  a  study  comparing  the  two  approaches  has  been  reported 
showing  each  of  these  to  have  its  area  of  superiority  [3].  Using  as 
improvement  factor  the  ratio 

F  -  (SIR)  /(SIR) 
o  e 

where  (SIR)o  is  the  output  signal-to-interference  ratio  at  the  array 
output  and  where  (SIR)^  is  Che  same  ratio  at  each  element  output,  we  find 
that  under  ideal  conditions — the  directionally  constrained  array  is 
accurately  pointed  and  the  LMS  canceler  array  is  equipped  with  a  perfect 
signal  reference  generator —  the  former  is  superior  when 


where  9^  is  interference  arrival  angle  relative  to  that  of  the  desired 
signal.  It  is  also  superior  when  the  inequality  is  reversed  if,  also, 
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the  interf erence-to-threshold  ratio  (ITR)  is  below  some  value  which 


decreases  as  ^  decreases.  (The  threshold  power  approximately  equals  the 
2 

noise  power  o^.)  Both  cancellers  suffer  under  conditions  of  inaccurate 
constraints;  that  is,  when  pointing  is  in  error  in  one  and  when  the 
reference  is  imperfect  in  the  other. 

An  alternative  which  suggests  itself  to  overcome  the  effects  of 
imperfections  of  these  kinds  is  a  hybrid  utilizing  both  types  of 
constraints.  One  can  be  expected  to  compensate  for  the  inadequacy  of 
the  other.  To  this  end  we  propose  the  scheme  shown  in  Figure  1.  The 
figure  shows  a  two-element  array  and  the  results  below  are  for  this 
case;  we  see  no  barrier  to  extending  the  scheme  to  multiple  elements. 

The  two  loops  in  the  figure,  opearting  on  the  two  array  ports  Z  and  A, 
function  as  a  least  mean  square  (IHS)  system  with  a  reference  signal 
extracted  from  the  array  output.  $  is  a  beam  steering  weight  which,  when 
properly  set,  generates  phase  coherent  signals  at  point  1  and  2  for  an 
EH  wave  arriving  from  a  prescribed  direction  representing  the  desired 
source.  Therefore,  $  sets  the  directional  constraint  and  Figure  1 
becomes  a  combinea  canceler. 

The  improvement  in  SIR  of  this  combined  canceler  can  be  shown 

to  be 

cos2 (9  /2)  |  l+2g  (1-F  )  1 1 1  2 sin6  (t'an(9  /2)-tan9  /2) ) 
p  _  _ I _ i _ ± _ i _ £ _  (i) 

cos2(Si/2)|  l+2g(l-F)s!s;2sin9g(tan(es/2)-tan(6i/2))  2 

where  F  and  F.  are  linear  oDerators  that  represent  the  effect  of  the 
si 

reference  loop  filter  on  the  desired  signal  and  on  the  interference 

2  ' 

spectral  contents,  respectively.  1/g  ■  (l+2c,.)/K  :  Pth.  Pth  stands  for 
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che  system's  power  threshold,  and  K  is  the  main  loop  d.c.  gain. 

A  w  d  ^  w^d 

9S  -  ®s]_"es'rM  ~ sia^s  and  "  9n“®i2  “  ~c~  sin’f|£»  where  d  is  the 

distance  between  the  array  elements,  and  u>^  are  the  radian  frequencies 

of  the  desired  and  interference  signals,  respectively,  c  is  che  speed 

of  electromagnetic  wave  and  i Ji  and  are  the  directions  of  arrival  of 

these  signals  with  respect  to  the  array  broadside. 

If  the  system  is  in  an  acquisition  mode,  the  reference  signal  loop 

generates  no  output  which  is  in  correlation  with  the  desired  signal.  In 

this  case  F  »  Fg  ■  0  (see  [5])  and  (1)  becomes  exactly  the  improvement 

factor  of  the  directionally  constrained  interference  canceler.  With  0 

at  a  value  close  to  that  needed  for  accurately  constraining  the  array 

Coward  che  desired  signal  we  can  expect  the  SIR  to  be  good  enough  to 

result  in  code  synchronization,  turning  the  system  to  a  code  tracking  mode. 

(The  code  acquisition  circuit  is  noc  shown  in  the  diagram.)  As  a  result, 

je 

J  sr 

Fj  remains  approximately  zero,  but  Fg  becomes  close  to  one  (F,.  -  ye 
where  y  and  9sr  are  the  amplitude  and  phase  of  the  filter’s  response). 

For  che  ideally  filtered  case;  Fg  ■  1,  (1)  then  becomes 

-  l+2g| I1 2sin6  (tan(9 . /2)-tan(9  /2)))2 

F  -  cos  (6  /2)[ - 5 - 1 - 5 - 1  (2) 

cos  (9i/2)  U; 

■  F  cos2(5  /2) 
o  s 

We  have  found  in  other  studies  [4  ]  that  the  sensitivity  to  an 
inaccurate  spatial  constraint  may  result  in  loss  of  SIR  in  the  order  of 
3  dB  for  0.1  degree  pointing  error.  This  result  shows,  however,  that 
with  perfect  reconstruction  of  the  reference  signal  the  imperfect  refer¬ 


ence  mainly  affects  the  improvement  in  che  factor  cos  (?s/2),  being 


1 


the  pointing  error.  This  effect  is  virtually  trivial  compared  to  that 
obtained  in  a  conventional  spatially  constrained  array. 

Another  potentially  useful  step  is  the  direct  correction  of  the 
pointing  error  via  the  processor  shown  in  Figure  2.  This  processor  contains 
the  following  control  loops: 

1.  The  main  weight  control  loop, 

2.  The  reference  signal  loop, 

3.  The  beam  steering  weight  control  loop 

4.  Reference  loop  phase  compensation 

The  first  and  second  loops  were  included  in  the  hybrid  scheme  of 
Figure  1.  The  hybrid  which  generates  sum  (I)  and  difference  (A)  signals 
eliminates  the  desired  signal  at  the  A  terminal  if,  indeed,  pointing  is 
accurate.  With  pointing  error,  this  will  not  be  the  case. 

The  purpose  of  the  third  loop  (the  beam  steering  weight  control  loop) 
is  to  adaptively  control  the  weight  $,  so  that  the  desired  signal  residue 
at  the  output  A  is  minimized.  This  is  done  by  correlating  (in  correlator 
#2)  the  output  of  the  reference  signal  loop  with  the  output  A,  of  the  hybrid. 

In  tracking,  the  output  of  the  reference  signal  loop  is  a  sufficiently 
clean  replica  of  the  desired  signal.  Ideally,  the  output  of  correlator 
number  (2)  will  settle  on  the  zero  value  (due  to  the  succeeding  integrator) 
only  when  the  output  A,  contains  no  power  that  is  correlated  with  the  desired 
signal.  In  a  practical  sense,  the  desired  signal  power  at  the  output  of 
hybrid  A  is  not  totally  eliminated,  but  rather  minimized  by  the  beam  steering 
loop  which  acts  in  the  IMS  fashion. 

In  the  acquisition  mode  the  output  of  the  reference  signal  is  very 
small  (due  to  the  filtering  of  both  the  spread  desired  signal  and  the 
interference).  To  prevent  output  of  correlator  number  2  from  incorrectly 
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affecting  the  steering  weight,  a  code  tracking  inhibit  circuit  is  used. 
That  is,  whenever  the  system  is  not  in  code  synchronism,  the  input  to 
integrator  number  2  is  inhibited,  disabling  the  automatic  beam  steering. 

In  such  a  case  (acquisition  mode  with  very  high  interference-to-signal 
ratio)  we  have  shown  (see  [3])  that  the  constrained  interference  canceler 
performs  as  a  power-inversion  device,  resulting  in  a  high  SIR,  and  leading 
back,  to  code  synchronization. 

The  fourth  loop  corrects  for  any  phase  shift  and  attenuation  caused 
by  the  BPF  [6). 
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The  technique  of  beamforming  on  a  beacon  to  organize  an  array  of 
arbitrarily  placed  sensors  is  described  in  the  literature  [1].  Essentially 
these  methods  involve  maximizing  the  array  output  to  the  signal  by  observ¬ 
ing  the  phase  of  the  arrivals  at  the  various  elements,  and  phase  shifting 
to  coherently  combine  them.  It  has  also  been  shown  that  once  the  array  is 
focused,  the  beam  can  be  steered  toward  another  source  by  open  loop  phase 
correction  using  the  steering  angle  and  nominal  array  element  position 
information.  The  angle  through  which  it  may  be  steered  with  a  loss  of 
response  of  less  than  1  dB  is 


where  X  is  the  received  signal  wavelength  and  0  is  the  rms  error  in  the 
knowledge  of  the  element  positions.  Typically,  at  X  -  1  cm,  a  1  cm 
uncertainty  in  element  position  will  allow 

A9  •  0.08  radian 

steering,  which  translates  to  about  1800  miles  on  the  ground. 

We  propose  the  following  approach  here  where  the  initial  focusing 
may  have  to  be  carried  out  in  the  presence  of  high  amplitude  interferers 
arbitrarily  located.  See  Figure  1.  The  beacon  will  be  assumed  adequately 
separated  from  potential  Interferers  so  that  they  will  not  influence  the  beam 
toward  the  beacon  when  it  is  formed.  The  beacon  is  assumed  to  be  of  strength 


o 


ARRAY  IN  SPACE 


BEACON  USER 

FIGURE  1.  BEAMFORMING  AND  SCANNING  IN  THE  PRESENCE  OF  INTERFERERS 

comparable  to  that  of  the  lnterferers  and  nodulated  In  a  manner  recognizable 
by  the  satellite  array.  We  will  assume  that  the  beacon’s  signal  Is  received 
separately  at  the  satellite  perhaps  via  a  conventional,  relatlvelv  low 
gain,  mechanically  steered  antenna.  (This  link  may  also  be  used  to  convey 
information  to  the  array  concerning  spread  spectrum  codes,  timing,  control 
signals,  etc.)  We  Imagine  the  array  deployed  on  a  large  ring  or  on  crossed 
ami  with  the  antenna  centered.  The  signal  so  received  will  be  used  as 
reference  input  to  the  array  operated  Initially  as  a  Widrow  LMS  processing 

* 

array.  The  steady  state  weights  generated  in  this  mode  are  given  closely  by 

*This  result  is  based  on  an  LMS  loop  with  large  open  loop  gain.  The  weights 

are  more  accurately  given  by  W  -  g(I+gM)-1  V*e_  where  g  is  the  open-loop  gain 
and  I  Is  the  unit  matrix.  r 


where  W  is  the  weight  vector 

M  *  is  the  inverse  of  the  covariance  matrix  of  array  element  outputs 
V  is  the  vector  of  signal  inputs 
e^  is  the  reference  signal 

The  bar  represents  time  averaging  over  an  interval  of  duration  T.  The  signal 
vector  is  comprised  of  components. 

v±(t)  -  cr(t)e^  S’  +  oq±(t) 

where  a(t)  is  the  signal  modulation  on  the  desired  carrier  arriving  from  the 
beacon 

$  is  the  carrier  phase  relative  to  some  reference 

is  a  complex  unit  amplitude  phasor  depending  on  the  arrival 
direction  of  the  signal  and  the  element  position 

o  . (t)  is  noise  and  interference 
The  reference  signal  will  have  the  form 
er  ■  ao(t)e^ 

where  a  is  an  arbitrary  real  constant.  Thus  the  vector  V*e^  will  have  components 

v^er  -  y  |  a|°2(t)|s1*dt  +  y  j  (t)e“^  0ql(Odt 

Assuming  a(t)  uncorrelated  with  the  noise  and  interference,  the  second  integral 

will  be  zero  if  T  is  made  large  enough  and  the  first  integral  will  be 
** 

a  P  Sj  where  P  is  a  measure  of  the  power  of  the  beacon  signal.  Thus 

-1  ’* 

W  «  a  P  M  S 

»* 

where  S  is  the  steering  vector  to  the  beacon.  If  the  receiver  simultaneously 


derives  the  covariance  matrix  M  then  forming  the  matrix  product  MW  one  has 

»* 

MW  -  a  P  S 

and  the  steering  vector  is  isolated.  Having  the  steering  vector  to  the  beacon, 
the  altered  vector,  S^,  which  will  point  the  array  toward  the  desired  source  is 
obtained  as  described  in  [1],  The  entire  system  is  represented  by  the  block 
diagram,  Figure  2. 


SIGNAL 

OUTPUT 


FIGURE  2.  BLOCK  DIAGRAM  OF  SELF-ORGANIZING  SYSTEM 
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To  summarize,  Che  array  initially  functions  as  an  adaptive  spatial 


filter  selecting  a  beacon  signal  from  a  multitude  of  signal  arrivals. 

The  steering  vector  inherent  in  the  adaptively-formed  element  weights  is 
isolated  and  altered  for  pointing  the  array  toward  a  desired  ground  user. 
Pointing  information  is  thus  generated  even  in  the  presence  of  unwanted 
signal  arrivals.  Questions  concerning  practical  implementation  and  the 
error  in  determining  the  pointing  vector  in  this  manner  remain  to  be 
dealt  with. 
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